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Abstract 


The  theory  of  Massera  and  Schaffer  relating  the  existence 
of  unique  almost  periodic  solutions  of  an  inhomogeneous  linear 
equation  to  an  exponential  dichotomy  for  the  homogeneous 
equation  has  been  completely  extended  to  discretizations  by 
a strongly  stable  difference  scheme.  In  addition  it  has  been 
shown  that  the  almost  periodic  sequence  solution  will  converge 
to  the  differential  equation  solution  at  a rate  O(k^)  where  p 
is  the  accuracy  of  the  scheme,  uniformly  in  t,  if  the  coeffi- 
cients are  sufficiently  smooth. 

The  preceding  theory  has  also  been  applied  to  a class 
of  exponentially  stable  partial  differential  equations  to 
which  one  can  apply  the  Hille-Yoshida  Theorem.  It  is  possible 
to  prove  the  existence  of  unique  almost  periodic  solutions  of 
the  inhomogeneous  equation  which  can  be  approximated  by 
almost  periodic  sequences  which  are  the  solutions  to  appropri- 
ate discretizations.  Two  methods  of  discretizations  are 
discussed;  the  strongly  stable  scheme  described  above  and 
the  Lax-Wendroff  scheme. 
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Introduction 


This  work  extends  certain  facets  of  the  theory  of 
Massera  and  Schaffer  [1]  (in  the  future  referred  to  as  M&S) 
relating  properties  of  the  solution  of  an  inhomogeneous 
ordinary  differential  equation  (ODE)  to  properties  of  the 
solutions  of  the  homogeneous  system.  In  Part  I we  extend 
this  theory  to  difference  approximations  of  the  ODE.  In 
Part  II  we  use  the  ODE  theory,  exclusively  in  the  stable 
case,  together  with  the  Hille-Yoshida  Theorem  to  obtain 
results  for  partial  differential  equations. 

The  study  of  admissibility  theory  and  dichotomy  theory 
for  difference  equations  was  first  done  by  Coffman  and 
Schaffer  [2],  This  work  differs  from  their  work  in  that  we 
are  concerned  with  the  preservation  of  admissibility  and 
dichotomy  properties  by  a difference  equation  which  is  used 
to  approximate  an  ODE  system  which  has  certain  of  these 
properties . 

First  let  us  review  the  theory  of  M&S.  Consider  the 
equations 

(0.1)  y = A ( t ) y 

(0. 2)  y = A(t)y  + f (t) 

Here  the  independent  variable  t ranges  over  the  whole  real 
line  and  for  each  t the  vector  y(t)  lies  in  a Banach 
space  E (which  may  be  infinite  dimensional).  A(t),  for 
each  fixed  t,  belongs  to  the  space  of  bounded  operators  on  E 
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which  we  will  denote  by  L(E,E).  Let  B and  D be  two 
Banach  spaces  of  a function  from  the  real  line  into  E. 

We  say  the  pair  (B,D)  is  admissible  for  (0.2)  if 

a)  V f e B 3 a unique  solution  y (to  (0.*2)  ) e D,  and 

(0.3) 

b)  fl  yll  D £ Kll  fll  B . 

This  definition  is  more  restrictive  than  that  of  MSS.  We  are 
not  going  to  deal  with  all  the  subtleties  of  their  theory, 
but  just  with  certain  important  parts.  It  is  pointed  out, 
however,  that  (0.3b)  is  actually  superfluous  (see  M&S, 

Chapter  5)  but  we  include  it  in  the  definition  of  admissi- 
bility for  simplicity. 

Observe  that  (0.3b)  states  that  if  we  write  the  y 
given  in  (0.3a)  as 

(0.4)  y = c(f) 

then  C,  which  is  obviously  a linear  operator  from  B into  D, 
is  also  bounded.  We  will  be  concerned  almost  exclusively 
with  the  case  that  D = B = A(E)  where  A is  the  space  of 
almost  periodic  (AP)  functions  with  range  in  E.  The  range 
in  this  notation  will  generally  be  omitted  and  we  will  simply 
write  A,  etc.,  if  no  confusion  might  arise.  At  times  we 
will  also  be  concerned  with  the  cases  B = D = L or 

OO 

6=0=0^  which  we  define  as  the  space  of  bounded  continu- 
ous functions  with  range  in  E.  Observe  that  A and 
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are  closed  subspaces  of  and  we  have  the  inclusions 


(0.5)  A C C C L 

00  00 

A(t)  will  always  be  an  AP  or  Lro  operator  function 
unless  stated  otherwise. 

Associated  with  the  concept  of  admissibility  we  define 
the  concept  of  an  exponential  dichotomy  for  (0.1).  Specifi- 
cally (0.1)  has  an  exponential  dichotomy  if  there  exist 
projections  P , P2=  I - P such  that  if  Y(t)  is  the 
fundamental  operator  solution  to  (0.1)  (Y(0)  = I,  see  M&S, 

Chapter  3) , the  following  estimates  hold  for  some  a > 0. 

a)  II  Y(t)P1Y_1(s)  II  < Ke“a(t“s)  , t <_  s 

(0.6) 

b)  II  Y (t)  P2y”  1 ( s ) II  < Ke~a(s_t)  , s >_  t . 

Note  that  this  is  not  the  general  definition  of  an 
exponential  dichotomy  given  in  M&S  (Chapter  8)  but  it  is 
equivalent  to  their  definition  when  E is  finite  dimensional. 
In  the  case  that  E is  infinite  dimensional  we  will  in 
general  only  be  concerned  with  the  stable  case;  that  is, 

P 1 = X'  P2  * °* 

Observe  that  P^  is  merely  the  projection  on  E+  , the 
space  of  initial  data  of  solutions  to  (0.1)  which  are  bounded 
for  t € [0,<»).  To  see  this  note  that  if  yg  were  in  the 
range  of  P2  and  if  the  corresponding  solution  y(s)  were 
bounded  for  s ^ 0 then,  by  setting  t = 0 in  (0.6b),  we 
obtain 
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(0.7) 


yQ  = Y 1(s)y(s)  = P2Y  1(s)y(s) 

and  yQ  = 0 follows  from  letting  s °°.  Similarly  P2  is 
just  the  projection  on  E_  , the  space  of  initial  data  which 
are  bounded  for  t £ (-oo,0].  An  exponential  dichotomy  merely 
states  that  E = E+  © E_  and  that  the  solutions,  in  addition 
to  being  bounded,  decay  exponentially. 

Now  it  is  a fundamental  result  of  M&S  that  if  A(t)  £ A 
then  (A, A)  is  admissible  for  (0.2)  iff  (0.1)  has  an  expon- 

ential dichotomy.  The  preceding  statement  is  also  true  if 
A is  replaced  everywhere  by  L^.  We  prove  here  the  easy  part 
of  the  result;  namely  if  we  have  an  exponential  dichotomy, 
then  we  have  admissibility.  If  f e Lot  define 


00 


(0.8) 

y ( t)  = c(f)  = | G(t,s) 

f(s)  ds 

where 

— 00 

G(t,s)  = YftjPjY-1 

(s)  , t > s 

(0.9) 

, -1 

= -Y(t)P2Y 

(s)  , s t . 

Observe  that  the  estimate 

(0.10)  I G ( t , s ) I - Ke_0ll  t_Sl 

holds  by  the  definition  of  an  exponential  dichotomy. 

Using  (0.10)  we  see  that  the  integral  in  (0.8)  exists, 
that  y satisfies  (0.2),,  and  that  y is  the  unique  Lot 
solution  to  (0.2)  and  in  fact 
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(0.11)  lyll  < — I fll  . 

■*  co  — (X  oo 

So  that  C,  defined  in  (0.8),  is  a bounded  operator  from 
Loo(E)  -+  ( E ) - To  see  that  C:  A -*■  A we  merely  let  T 

be  a common  e-almost  period  for  f and  A.  Then  w(t)  = y(t+T) 

- y(t)  is  the  unique  solution  to  (0.2)  with 

inhomogeneous  term 

(0.12)  [A(t+T)-A(t) ]y (t+T)  + f(t+T)  - f(t)  =0(e) 

We  then  use  (0.11)  to  obtain 
(0.13)  II  (y ( t+T)  - y(t)  )Wm  = 0(e)  . 

The  proof  of  the  converse  is  more  difficult  and  can  be 
found  in  M&S,  Chapter  10.  A simplified  proof  valid  only  in 
the  finite  dimensional  case  is  given  in  the  Appendix. 

In  Section  1 we  will  define  almost  periodic  sequences, 
which  will  be  the  type  of  solution  we  will  be  searching  for. 
Section  2 is  the  most  important  of  this  work.  Here,  after 
discussing  the  properties  of  the  strongly  stable  difference 
schemes  we  will  be  using,  we  will  introduce  a transformation 
which  will  separate  out  the  roots  of  the  scheme  inside  the 
unit  circle  and  permit  us  to  work  with  the  Euler  1-step  scheme. 
This  technique  was  originally  developed  by  Engquist  [3] 
although  the  author  was  not  aware  of  his  work  when  the 
formulation  given  in  Section  2 was  developed.  Engquist' s 
results  will  be  discussed  more  thoroughly  at  the  end  of 
Section  3. 


5 


i ’ 

li 


m>r. 
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In  Section  3 we  settle  the  question  of  admissibility 
when  the  homogeneous  system  is  exponentially  stable.  Using 
an  inequality  which  is  an  exact  discrete  analogue  to  the 
Gronwall  inequality,  we  will  show  that  the  homogeneous 
difference  equation  is  also  exponentially  stable.  From 
there  it  will  be  a simple  matter  to  obtain,  for  any  suffi- 
ciently small  time  step  k,  the  existence  of  an  almost 
periodic  sequence  as  a solution  to  the  inhomogeneous 
difference  equation.  Furthermore  the  sequences  converge 
uniformly  to  the  unique  AP  solution  to  the  ODE  with  a 
uniform  error  O(k^)  where  p is  the  order  of  accuracy  of 
the  scheme.  We  will  also  show  that  this  solution  can  in 
fact  be  calculated,  i.e.  it  is  stable  under  roundoff 
errors  and  errors  in  initial  data. 

In  Sections  4 and  5 we  deal  with  the  case  that  the 
homogeneous  system  has  a general  exponential  dichotomy. 

We  will  show  that  the  corresponding  inhomogeneous  difference 
equation  also  has  an  exponential  dichotomy.  This  is  only 
of  theoretical  interest  as  the  solution  will  no  longer  be 
stable  under  roundoff  errors  or  errors  in  initial  data. 

In  Section  6 we  will  deal  with  some  miscellaneous  topics, 
especially  the  convergence  of  the  mean  value  of  the  AP 
sequence  to  the  mean  value  of  the  AP  solution,  and  also  the 
weakly  nonlinear  case. 

In  Part  II  we  extend  this  theory  to  a simple  class  of 
partial  differential  equations  which  can  be  written  as  an 
evolution  equation 
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(0.14) 


y 


[B— 6 ] y + f 


l 

f 


where  B is  an  unbounded  operator  which  satisfies  the 
conditions  of  the  Hille-Yoshida  Theorem  and  6 (t)  is  an 
AP  function  such  that  the  homogeneous  system  is  exponentially 
stable.  We  can  obtain  a unique  AP  solution  to  (0.14)  by 
using  the  formula  which  would  be  valid  if  B were  bounded 
and  then  showing  that  under  mild  restrictions  on  f the 
resultant  function  does  in  fact  satisfy  (0.14). 

In  Section  8 we  introduce  a family  of  bounded  operators 
B^  which  are  spatial  discretizations  to  B.  We  construct 
functions  yh  which  are  the  unique  AP  solutions  to 

(0.15)  yh  = [Bh~6]yh  + f . 

We  will  give  conditions  to  insure  that 

(0.16)  II  (y-yh)  I oo=  0(hj)  , 

where  j is  the  order  of  the  approximation  of  B^  to  B. 

In  Section  9 we  apply  the  theory  of  Part  I to  (0.15)  to 
obtain  an  AP  sequence  yn.^  which  approximates  y^. 

We  have  however  the  unfortunate  restriction 


(0.17) 

where  II  B,  II 


2m 


= 0(1) 


m. 


0(l/h  ).  A more  favorable  result 


(0.18)  = 0(1) 

h 

is  obtained  in  Section  10  when  using  the  Lax-Wendroff  scheme. 
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In  Section  11  we  will  consider  the  extension  of  these 
results  in  the  case  that  the  operator  B is  perturbed  by 
some  bounded  AP  perturbation  D(t). 

Finally  in  the  Appendix  we  will  give  a proof  of  the 
basic  ODE  theorem  that  admissibility  is  equivalent  to  the 
existence  of  an  exponential  dichotomy.  The  proof  is  valid 
only  in  the  finite  dimensional  case  but  is  simpler  than  the 
proof  given  in  Mf.S  and  is  also  simpler  than  a finite 
dimensional  proof  to  be  found  in  Coppel  [4]. 
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1 • Almost  Periodic  Sequences 


Our  first  task  is  to  introduce  the  discrete  analogue 
of  an  AP  function.  Following  Corduneanu  [5]  we  define  an 
almost  periodic  (AP)  sequence  an  with  range  in  E as 
follows : 

is  AP  iff  given  e > 0 there  exists  a length  L(e) 
(a  positive  integer)  such  that  in  any  sequence  of  L 
consecutive  integers  there  exists  an  N such  that 


(1.1) 


a II  < 
n 00 


e 


where  the  sup  in  (1.1)  is  taken  over  n.  As  shown  in 
Corduneanu  (page  45)  this  is  equivalent  to  normality  i.e. 
given  any  sequence  of  integers  the  sequence 

bn;i  = an+N.  wil1  have  a uniformly  convergent  subsequence. 
Although  the  proof  given  in  Corduneanu  is  stated  for  scalar 
valued  sequences  this  proof  is  obviously  valid  if  the  range 
is  any  Banach  space  E. 

If  we  define  the  space  L^(E)  as  the  Banach  space  of 

bounded  sequences  with  range  in  E then  the  AP  sequences  An(E) 

form  a closed  subspace  of  L^(E).  As  usual  the  argument  E 

will  be  omitted  when  no  confusion  can  arise. 

We  point  out  that  normality  can  be  used,  exactly  as  in 

the  continuous  case,  to  show  that  for  any  finite  set  of  AP 
1 r 

sequences  a , . . . , an  with  range  in  possibly  different  spaces 
1 r 

E ,...,E  and  for  any  e > 0,  there  is  always  a length  L(e) 


such  that  in  any  interval  of  length  L we  can  find  a common 
e-almost  period. 

Finally  we  note  that  if  f(t)  is  an  AP  function  then 
the  sequence  f = f(nk)  is  an  AP  sequence  for  any  real  k. 
The  converse  is  also  true  as  shown  by  Corduneanu  (page  47) 
but  we  shall  not  use  that. 


2 . Properties  of  the  Difference  Scheme 


We  consider  linear  £-step,  strongly  stable  schemes 
described  as  follows. 


• AJ 

(2.1)  l a .y  . = k 7 3 .y  . = k 7 3 . [A  ,.y  ,.+f  .] 


j=0 


j=0 


D n+] 


j = 0 


j n+3Jn+]  n4] 


Here  we  assume  that  we  are  discretizing  the  ODE 


(2.2) 


y = A (t)  y + f(t) 


and  yn  = y(nk),  An  = A(nk),  f = f(nk)  where  k is  the 
time  step. 

Associated  with  (2.1)  we  have  the  polynomials 
l l 

(2.3)  p(x)  = l a.x-1  , a(x)  = £ 3.x-1  . 

j=0  j =0  J 

It  is  well  known  (see  Dahlquist  [6]  or  Heinrici  [7])  that 
consistency  implies  that  x = 1 is  a simple  root  of  p(x)  =0 
and  that 

(2.4)  p'  (1)  = a(l)  = 1 

where  we  have  normalized  the  coefficients  so  that  the  common 
value  in  (2.4)  is  one. 

Strong  stability  of  the  scheme  is  achieved  by  restricting 
the  size  of  the  other  £-1  roots  of  p(x)  = 0.  Specifically, 
if  we  number  these  roots  xu  , u = then  we  require 
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that  there  be  a positive  number  0 < 1,  such  that 


(2*5)  = 1 , jxj  < 8 < 1 , u = 2 . 

Finally,  we  make  two  further  assumptions. 

(2*6)  xu  distinct  , u = 2,...,£  , 

(2.7)  xu  i 0 , u = 2, . . . , SL  . 

/ 

(2.7)  is  necessary  because  to  get  an  AP  solution,  the 
difference  equation  should  be  solved  backwards  and  forwards. 
This  assumption  can  be  removed  in  the  important  case  that 
the  homogeneous  ODE  system  is  exponentially  stable.  Condition 
(2.6)  can  be  removed  in  all  cases  and  is  included  here  only 
to  simplify  the  following  proofs.  The  removal  of  these 
conditions  will  be  discussed  in  Section  6. 

Now  in  working  with  a multistep  scheme,  the  standard 
procedure  is  to  convert  it  into  a one-step  scheme.  To  do 
this  we  define  the  space  = E><Ex...xE.  (We  will  usually 
write  vectors  in  E^  in  column  vector  form. ) We  give  E^  the 
norm  inherited  from  this  definition,  namely  if  w e E„  and 

X/ 

r 1 

y 

(2.8)  w = 

• 

a 

y 

then 

(2.9)  II  wll  ^ = max  Hy^i  . 

. hj 

l 

Here  we  have  explicitly  indicated  the  E norm  in  (2.9). 
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We  now  consider  the  discretization  of  the  linear 


inhomogeneous  system  (2.2).  If  wn  = 


(2.10) 


Yn+£-l 


1 ^n 


then  we  get 


w .=Cw  + kf 
n+1  n n n 


Here  if  we  define  z = k A , 

n n 


(2.11)  cn  C(zn, . . . ,zn+£)  - 


(a£I~B£Zn+£)  (eA-lZn+A-l“Ct£-lI) 


0 

I 


■•(V“B£Zn+A,"1(B0Vo0I) 

0 


. . 0 


0 10 


f = 
n 


‘“t1 


-1 


£ 

l 3 

j=0 


j fn+j 


Observe  that  Cn  € L(E£,E£)  and  will  be  just  an  m£xm£ 
matrix  in  the  case  that  E is  an  m-dimensional  space  (in  which 
case  it  is  called  the  companion  matrix) . Note  also  that  since 
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A is  uniformly  bounded  in  t,  the  implicit  term  (a.l-g  z ) 

£ £ n+£ 


can  be  inverted,  for  k sufficiently  small,  uniformly  in  n, 


We  finally  point  out  that  if  A and  f are  AP  (L  ) then  C 

00  r 


311(1  \ are  AP  (l£)  and  also  that  C given  in  (2.11)  is 


smooth  function  (in  the  Frechet  sense)  of  its  £+1  arguments. 


Now  with  A a constant  and  z = kA  e L(E,E)  we  consider 


for  small  II  zll  the  homogeneous  difference  scheme 


<2'12)  wn+l  = U(2)wn 


U(z)  = (a£I-Btz)‘1(Bl_1z-a1,_1I)  ...  (a£I-B£z) _1  (B0z-a0I ) 


U(z)  is  a mapping,  defined  for  small  II  zll  , from  L(E,E)  into 


L^E£'E£^’  Observe  the  following  properties  of  U(z). 


(2.13)  (a)  U ( z ) is  a smooth  (in  the  Frechet  sense)  function  of 


(b)  U (0)  has  eigenvalues  exactly  with  eigenspaces 


Eu  , where 


x y 
u 


x y 
u 1 


y e e 


(c)  E = © J E . 

u=l  u 


P 


Now  (a)  and  (b)  should  be  clear.  To  see  (c)  , note  that  the 


spaces  Eu  are  closed  and  also  that  any  two  have  only  zero 


/I 

m common.  Let  w = .*  and  suppose  w has  an  expansion 

X£ 


(2.14) 


£ 

w = l 

u=l 


r £-1  ^ 

'x  y 
u . u 


x y 
u u 


v y 


u 


Then  if  we  define  the  vector  w to  be 
(2. 15)  w = V w 


r Yll 

• we  can  write 

y£J 


V = x 


£-1 


£-1 

c£ 


X1  ‘ ‘ ' X£ 

1 ...  1 


i.e.  V is  just  the  Vandermonde  matrix  associated  with  the 
distinct  numbers  x^,...,x^.  (2.13c)  now  follows  directly 

from  the  invertibility  of  V. 

The  representations  (2.14)  and  (2.15)  enable  us  to 
define  an  equivalent  norm  on  the  space  . Specifically  if 
w is  expressed  as  in  (2.14)  we  define 


(2.16)  liwnu  = max  |j  n 

u u E 

That  II  II  is  equivalent  to  II  II  ^ follows  immediately  from 
the  representation  (2.15)  and  the  invertibility  of  V.  Note 
chat  the  equivalence  of  these  norms  implies  the  equivalence 
of  the  operator  norms  they  induce  on  L(E£,E£).  Of  course 
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this  paragraph  is  superfluous  in  the  case  that  E is  finite 
dimensional. 

For  future  use  we  point  out  that  (2.13b,c)  imply 

that  U l(z)  exists  for  small  II  zll  and  is  smooth  in  z. 

Now  define  the  space  E^  = © £ E . We  have  E = E.8E1 

u>l  u 1 1 

and  we  note  that  E^  is  canonically  isomorphic  to  E.  At  times 
we  will  identify  E^  with  E but  this  should  not  cause  any 
confusion.  With  respect  to  the  decomposition  E = E^©  E^  we 
see  that  U(0)  is  in  block  diagonal  form;  symbolically, 

rD (0)  0 

(2.17)  U (0)  = 

1 0 B ( 0) 

here  D(0)  : E^  E1  and  is  the  identity,  while  B(0):  E^E1 

and  II B ( 0 ) II  U < 0.  This  is  the  operator  norm  induced  on 
k (E £ ' E£ > by  the  11  11  U norm  on  E^  and  follows  from  the  fact 
that  B(o)  is  just  multiplication  by  xu  on  the  space  Eu« 

We  can  now  state  the  fundamental  theorem  of  this  section. 

Theorem  1.  For  small  II  zll  there  exists  an  operator 
T(z):  L(E,E)  - L(E£,E£)  such  that 

(a)  T(0)  = I (Identity  on  E^ ) 

(b)  T(z),  T ■*'(z)  are  smooth  in  z 

(c)  L(z)  = T ^(z)U(z)T(z)  is  in  block  diagonal  form 
with  respect  to  the  decomposition  E^  = E^  © E^. 

(d)  Writing  L(z)  symbolically  as 


then  B ( z)  | z=0  = B (0 ) (from  (2.17))  and  II  B ( z ) II  U < 0 while 
D(z)  has  an  expansion 

(2.18)  D ( z)  = I + z + 0(11  z II  2 ) 


where  is  identified  with  E. 


Theorem  1 is  basically  trivial  and  the  proof  involves 
familiar  arguments.  First  we  note  that  U(z)  is  smoothly 
invertible  for  small  I zll  . (The  restriction  "for  small  II  zll  " 
will  not  be  stated  explicitly  in  the  future.) 

Next  we  observe  that  (WI  - U ( z ) ) “ 1 exists  for  w in  a 
small  annulus  around  the  circle  |w|  = 0 and  this  holds 
uniformly  in  z.  This  follows  from  (2.13b)  and  the  geometric 
series . 

Now  define  the  projections 

(2.19)  P(z)  = ~ j j (wi-u(z))  1 dw  , Q ( z)  = I - P(z)  . 

| w | =0 

That  P and  Q are  projections  is  a familiar  result  which 
follows  from  the  resolvent  identity, 


(2.20) 


(w1I-U(z))  1-  (w2I-U(z))_] 


= (w2~w1)  (w1I-U(z)  ) 1(w2I-  U ( z)  ) 1 , 


calculating  P by  integrating  around  two  slightly  different 
circles,  and  interchanging  the  order  of  integration.  It 
follows  from  the  construction  that  P(z)  and  Q(z)  commute 
with  U(z)  and  that  they  are  smooth  functions  of  z. 
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It  is  clear  by  applying  P(0)  to  an  arbitrary  vector 
w G E£  311(1  usin9  the  expansion  given  in  (2.15)  that  P(0) 
is  exactly  the  projection  onto  E1  along  and  Q(0)  is  the 
projection  onto  E1  along  E1. 

Now  define 

(2.21)  T ( z)  = P (z) P (0)  + Q (z)Q (0)  . 

Clearly  T(0)  = I and  T(z)  is  smooth  whence  T_1(z)  exists 
and  is  also  smooth.  Note  that  this  holds  in  either  of  the 
norms  II  II  U or  II  II  £ on  L(E£,E£).  Let 

(2.22)  L ( z)  = T~1(z)U(z)T(z)  . 

We  claim  L is  in  block  diagonal  form  with  respect  to 
the  decomposition  E£  = E±  © E1  . This  follows  directly 
from  the  fact  that  P(z)  and  Q(z)  commute  with  U(z).  In 
fact  let  x G E^  ; then 

(2.23)  U (z)  T (z)  x = U ( z)  Q ( z)  x = Q(z)U(z)x  = w (say). 

Now  if  y = T 1 ( z)  w and  y = yx  + y1  with  y G e1 
and  y ^ G e ^ we  have 


j 


(2.24)  w = Q ( z)  U ( z)  x = T ( z)  y = Q(z)y1  + P(z)y1 
whence 

(2.25)  P(z)y^  = T(z)y1'  = 0 

and  so  y1  = 0 by  the  invertibility  of  T(z).  Thus  L(z) 
maps  E ^ ■>  E ^ and  a similar  argument  shows  L(z)  : E1  -*•  E1. 
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If  we  write  {using  the  notation  of  (2.17))  L(z)  as 


(2.26)  L ( z) 


D { z)  0 D ( z)  : ->  E ^ 

0 B (z)  B (z)  : E1  -*■  E1 


then  we  see  that  D(z)  and  B(z)  are  smooth  in  z;  D(0)  = I 
and  II  B ( 0 ) II  U < 6.  It  follows  immediately  by  continuity 
that 


(2.27) 


II  B ( z ) II  U < 0 


for  small  II  zll  . Now  (2.27)  will  imply  that  given  z^,...,z^ 
with  II  z .11  small  then 

(2.28)  II  FT  B(z.)ll£  < k ed 

i = l 1 

for  a constant  K independent  of  the  z i and  N. 

It  only  remains  to  study  the  first  order  structure  of 
D(z)  in  order  to  obtain  (2.18).  Expanding  D(z)  about  z = 0 
we  can  write 

(2.29)  D ( z)  = I + C ( z)  + 0(11  zll  2 ) 


where  C is  the  Frechet  derivative  of  D(z)  at  z = 0 and  is 
a bounded  linear  map  from  L(E,E)  into  L(E^,E^)  which  we 
identify  with  L(E,E).  It  is  necessary  to  show  that  C is 
the  identity,  and  as  one  might  expect  this  will  follow  very 
easily  from  consistency. 

Let  A be  an  arbitrary  element  in  L (E,E)  and  let  z = kA, 
a'  = C (A)  . For  y G E and  identifying  E^  with  E,  we 
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calculate  the  following  limit 


(2.30)  lim  D(kA)  y = (I  + kA  ' +0(k2))ny 
k+0 

nk=  t 

for  any  fixed  t > 0. 

As  one  would  expect  this  limit  is  etA  y because 
the  difference  scheme  expressed  in  (2.30)  is  consistent 
with  the  ODE  y = A'y. 

More  precisely  since 

(2.31)  « D ( k A ) II  < 1 + kR 


for  small  k,  where  R is  some  fixed  constant,  it  follows 
that 


(2. 32) 


D(kA)N«  < eRT  , 


Nk  < T,  0 < t < T . 


Let  xn  s°lve  the  difference  equation 


(2.33) 


xn+i  = D(kA)xn 


x0  = y 


i.e.  xn  D(kA)  y.  Now  if  x^  - e ny  then  xn  satisfies 
the  perturbed  difference  equation 


(2. 34) 


xn+i  = [D(kA)  + 0(k2)  ]xn 


where  the  term  0(k  ) is  uniform  for  nk  £ T.  Equation  (2.34) 
follows  directly  from  the  fact  that  x(t)  satisfies  the  ODE 


(2.35) 


x = a'  x . 
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fW*1' 


Letting  wn  - xn  ” xn  we  can  show  easily  enough  that 


(2. 36) 


lwnll  = 0 (k)  , 


nk  < T . 


In  fact  from  (2.33)  and  (2.34)  we  obtain 


(2.37) 


* = k l D(kA)n-]  0 ( 1)  w . 

j=l 


If  we  use  (2.32)  to  bound  the  powers  of  D(kA)  and  thus  II  x II  , 
we  then  obtain  (2.36)  and  letting  k 0 we  obtain 


(2.38) 


lim  D (kA)  y = e y 

n-*co 

nk=t 


Recalling  the  identification  of  E.  with  E we  define  the 
Y _ i 

vector  y = : whence  y e E.. 

y 1 


Thus  (2.38)  can  be  expressed 


(2.39) 


lim  L(kA)ny  = 
k+0 
nk=  t 


r A t 
e y 


A 1 1 
e y 


Now  we  are  going  to  show  that  A’  = A by  calculating  this 
limit  directly  from  the  definition  of  L (kA) . In  fact  since 
II  U ( 0 ) II  U =1  it  follows  that  II  U ( k A ) n II  u is  uniformly  bounded 
for  nk  £ T (and  of  course  this  is  also  true  if  we  replace 
II  II U by  II  11^).  Then  writing 

(2.40)  L(kA)n  y = T-1(kA)  U(kA)n  T (kA)  y , 


and  letting  k 0 such  that  nk  = t,  using  the  continuity 
of  T(z)  and  the  fact  that  T(0)  = I together  with  the 
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uniform  boundedness  of  the  powers  U(kA)nf 

(2.41)  lim  L(kA)n  y = lim  U(kA)n  y 

k-*0  k+0 

nk  = t nk=t 

the  last  equality  follows  from  the  fact  that  a consistent 
and  stable  difference  scheme  is  convergent. 

Since  A was  arbitrary,  (2.41)  together  with  (2.39)  shows 
that  A = a',  which  proves  that  the  map  C is  the  identity 
and  this  completes  the  proof  of  Theorem  1. 


we  are  left  with 

At 
e y 


At 
e y 
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:w  -'.H.  ----- 


... 


3.  The  Exponentially  Stable  Case 


We  now  assume  that  the  homogeneous  system  (0.1)  is 
exponentially  stable.  Thus,  if  Y(t)  is  the  fundamental 
operator  solution  to  (0.1)  (Y  (0 ) = I),  then  the  estimate 

(3.1)  II  Y(t)Y_1(s)H  < Ke"a(t_s)  , t > s , 


will  hold  with  positive  constants  K and  a.  In  the  sequel 
we  will  use  K as  a generic  positive  constant  so  it  will 
appear  in  contexts  other  than  (3.1). 

We  wish  to  prove  a similar  estimate  for  the  homogeneous 
difference  equation  (see  (2.10),  (2.11)) 


(3.2) 


w , . = C w 
n+1  n n 


Now  Cn  is  invertible  for  all  n by  (2.7).  Then  (3.2) 
has  a fundamental  solution  Wn  (Wg  = I)  which  is  simply 

n-1 

W = TT  C.  , n > 0 , 

i=0  1 

(3.  3) 

-n  _i 

W = TT  C.  , n > 0 . 

i=-l  1 


It  can  easily  be  verified  that  the  unique  solution 
to  (3.2),  w^  , given  initial  data  w^  is  simply 


(3.4) 


w = W W . ^w  . 
n n d j 


Note  that  for  n > j,  wnWj^  is  simply 


'.  r’5**  /w 


(3.5) 


W W . 
n 3 


n-1 


C. 

1 


n > j . 


In  the  stable  case  we  will  only  use 
and  we  can  see  that  the  invertibility 
in  this  case  and  assumption  (2.7)  is 
here . 


W W . 1 for  n > i 
n 1 - J 

of  is  not  required 

therefore  not  required 


We  intend  to  show  that  for  small  k (this  will  not  be 
stated  explicitly  in  the  future)  (3.1)  implies 

-a1k(n-j) 

(3.6)  IIWnWjll  <_  K e 1 , n > j , 


where  K is  used  as  a generic  constant  and  can  be  taken 
independently  of  k,  while  = a + 0(k).  Thus  since  a > 0, 
can  also  be  taken  independent  of  k for  small  k. 

To  attain  (3.6)  we  suppose  A(t)  is  C1.  All  derivative 

that  are  assumed  will  always  be  required  to  be  AP  or  L . This 

00 

will  not  be  stated  explicitly  in  the  future.  The  smoothness 
condition  on  A will  be  removed  in  Section  6. 

Expanding  An+j  about  A n we  can  write  (3.2)  as 


s 


(3.  7) 


w 


n+1 


= U (kA  ) w 
n n 


+ k 0 (1)  w 


n 


Here  the  0(1)  term  is  uniform  in  n and  the  operator  U is 

given  by  (2.12).  We  are  going  to  show  that  perturbations  of 
2 

the  type  k 0(1)  preserve  exponential  stability.  We  thus 
consider  the  unperturbed  system 


(3.8) 


wn+l  = U(kVWn 
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If  we  let  v = T(kA  )w  , where  the  operator  T is 
n n n 


introduced  from  Theorem  1,  and  note  that  by  the  smoothness 


of  T and  A we  have  T(kAn+1)  = T (kA^)  + k 0(1),  we  see  that 


we  should  consider  the  difference  equation 


(3.9) 


v . = L (kA  ) v 
n+1  n n 


where  the  operator  L has  been  studied  in  Theorem  1. 


If  denotes  the  fundamental  solution  to  (3.9)  then 


since  L is  block  diagonal  we  see  that  Vn  itself  is  also 


block  diagonal.  If  we  write  Vn  in  the  form 


(3.10) 


V = 
n 


0 V 


we  can  see  that  the  second  part,  the  contribution  from  the 


roots  inside  the  unit  circle,  causes  no  difficulty.  In  fact 


for  n > j we  have 


(3.11) 


? o-1  n~l 

V V . = 1 T B (kA.  ) 

n i \ \ l 

J i=D 


v2v?  1ii  < k en-^ 

n 3 


where  we  have  used  (2.28)  and  K as  a generic  constant. 


Since  9 is  a fixed  number  less  than  1 we  see  immediately 


that  (3.6)  holds  for  V where  a,  can  in  fact  be  taken 

n 1 


equal  to  a for  small  k. 


We  must  now  deal  with  Vn-  Neglecting  perturbations 


of  k 0(1)  and  identifying  the  space  with  E we 
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consider  the  difference  equation 


(3-12>  xn+i  = d + kAn)xn 

i.e.  we  have  reduced  the  problem  to  the  study  of  the  Euler 
1-step  scheme. 

We  will  compare  Xn  , the  fundamental  solution  to  (3.12), 
with  Yn  (=  Y(nk))  where  Y is  the  fundamental  solution  to 
the  homogeneous  equation  (0.1). 

Let  H(t,s)  = Y(t)Y  ^(s).  As  a function  of  t, 

H satisfies  the  equations 

H = A ( t)  H , 

(3*13)  H = A2(t)H  + A ( t)  H , 

H ( t , t ) = I 


H 


Integrating  (3.13)  we  obtain 


H ( (n+1) k ,nk)  = I + k H(nk+0k,nk)  d0 


(3.14) 


1 1 


= I + kA(nk)  + k‘ 


0 d0  d(J)  H (nk  + <j)0k,  nk) 


0 0 


Now  since  A is  bounded  it  is  a standard  result  (M&S, 
Theorem  3.1  C) , that 


(3.15)  II H ( t , s ) II  < K1(R)  for  |t-s|  <R  , 


where  R is  any  positive  number.  This  does  not  require 
exponential  stability  and  in  fact  is  an  immediate  consequence 
of  Gronwall's  inequality.  It  now  follows  from  the  second 
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equation  in  (3.13)  that  the  coefficient  of  k2  in  (3.14)  is 
bounded  uniformly  in  n and  k.  This  implies  that  is 

the  fundamental  solution  operator  to  an  equation 

(3.16)  Yn+1  = [I+kDn(k)]Yn  = [I+kAn  + k20(l)]Yn 

where  the  0(1)  term  is  uniform  in  n.  This  means  that  we 
can  regard  the  x^  equation  (3.12)  as  a perturbation  of  the 
Y^  equation.  We  point  out  that  the  reduction  we  have  obtained 
does  not  use  exponential  stability  and  so  is  valid  in  the  case 
that  (0.1)  has  a general  exponential  dichotomy.  If  we  use 
the  stability,  however,  we  see  that  Y satisfies  the 
estimate 

(3.17)  H YnY j 1,1  £K  e_ctk(n~j)  , n > j . 

We  must  now  prove  the  proposition  that  exponential  stability 

2 

is  preserved  under  k perturbations. 

Theorem  2.  Consider  two  difference  equations 

(3.18)  yn+1  = R(n,k)yn 

(3.19)  x , , = R(n,k)x  + kS(n,k)x 

n+1  n n 

defined  for  small  k.  Suppose  that  (3.18)  is  exponentially 
stable,  that  is  there  exists  constants  K and  a independent 
of  k such  that 

(3.20)  II  Y Y_1ll  < K e-ak  (n“  j)  n > j , 

n j — ' — J 
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where  Yn  is  the  fundamental  solution  to  (3.18).  We  point 

out  that  we  do  not  require  R(n,k)  to  be  invertible  since 

we  only  study  the  solution  yj’1  for  n > j.  In  that  case 

Yn  is  not  invertible  and  we  must  use  (3.5)  in  place  of  Y Y-1 

n j ' 

but  we  will  retain  this  notation  for  simplicity. 

Under  this  hypothesis  there  exists  an  eQ  such  that  if 
for  small  k. 


(3.21) 


S (n ,k) 


< £ < p 

- - Lo 


where  the  sup  in  (3.21)  is  over  n,  then  (3.19)  is  exponentially 
stable  and  if  Xn  is  the  fundamental  solution  to  (3.19)  (the 
dependence  on  k has  been  suppressed)  the  estimate 


(3.22) 


~ ~_i  -a.k(n-j) 

Vj  " < v 1 


where  K^.  can  be  taken  independent  of  k and  = a + 0(e). 

Before  proving  this  we  note  that  the  case  where  S(n,k) 

= 0(k)  is  automatically  covered.  The  more  general  formula- 
tion will  be  used  in  Section  6 and  more  importantly  in  Part  II. 
Note  also  that  it  is  a discrete  analogue  to  the  Gronwall 
inequality  and  the  proof  is  in  fact  immediately  suggested  by 
the  proof  of  the  Gronwall  inequality.  This  theorem  is  equivalent 
to  a lemma  of  Engquist  itself  based  on  a theorem  of  Strang  [8], 
but  the  proof  given  here  is  simpler  than  Strang's  proof  and 
much  more  suggestive  of  the  Gronwall  inequality. 

To  prove  Theorem  2 we  first  note  that  for  n > j we  have 


(3.23) 


Vj1= 


Y YT1  + k 
n j 


I 


£=j  + l 


Y yT1 
n 3 


SL-1 
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This  analogue  of  the  variation  of  constants  formula  can  be 
verified  immediately.  Using  (3.20)  we  get 


(3.24)  II  X X . 

n 3 


"In 


< K e-otk  (n-3 ) + kK  eak£  V ^-ak(n-Jl)  „ ~ —1 


l e 
£=j 


X . X . II  . 
I 3 


Here  K is  the  constant  of  (3.20).  If  we  redefine  K as 

K eak  (for  small  k)  and  define  v . = eakn  II  x xT  1l!  we 

n / 3 n 3 

can  write 


n-1 


(3.25)  v . < K + kKe  l v.  . E s . . , 
n ' 3 ~ A,3  n-1,]  ' 


n > j 


From  (3.25)  we  obtain 


(3.26) 

whence 


v 


n»3 


S .-s  , 

= n/3  n-1 , i 

kKe 


< s . . 

- n-1, 3 


(3.27)  sn^.  < (l+eKk)n-j  s.  j < eeKk(n-j)  , 


where  K is  defined  so  that  s . . < K.  independent  of  k 

and  j (from  (3.25)  specialized  to  n = j+1  and  the  fact  that 

v . . = 1)  . 

3 • 3 

Inequality  (3.27)  together  with  (3.25)  and  the  defini- 
tion of  vn  f j yields  (3.22)  immediately  and  thus  completes 
the  proof  of  Theorem  2 . 

Returning  to  our  specific  case,  it  should  be  clear  that 
after  several  applications  of  Theorem  2 we  can  prove 
exponential  stability  for  v defined  as  T(kAn)w  (see  (3.8) 


and  the  following  paragraph) . Since  in  terms  of  fundamental 
solution  operators  we  have 

(3. 28)  W W = T_1{kAn)V  V^TfkA.) 

n 1 n ] j' 

we  obtain  exponential  stability  for  (3.8)  and  another  applica- 
tion of  Theorem  2 yields  exponential  stability  for  the  full 
homogeneous  system  (3.2).  It  is  a simple  matter  to  go  from 
exponential  stability  to  admissibility  for  the  inhomogeneous 
difference  equation 


w = C w + kg 
n+1  n n yn 

We  first  define  admissibility  in  the  obvious  way. 

The  pair  (i£,l£)  will  be  admissible  for  (3.29),  for  small 
k,  iff  for  any  sequence  gn  in  l£  (3.29)  has  a unique  solu- 
tion wn  in  i£  and  this  assignment  is  a bounded  mapping, 
i.e.  there  exists  a K (independent  of  k)  such  that 


(3. 30) 


Vo  i K" g. 


n 


An  entirely  analogous  definition  holds  for  the  admissibility 
of  the  pair  (A  ,A  ) . As  in  the  ODE  case  we  will  show  (Ln,Ln) 
is  admissible  and  then  show  that  the  solution  is  AP  if  the 
coefficients  are. 

To  show  (L^l/J)  admissibility  we  merely  write  the  solutioi 


n 


w - k I W W . ^ g . , 
n j=-oo  n 3 3-1 


(3. 31) 

Exponential  stability  implies  the  convergence  of  this  series. 
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which  can  be  immediately  verified  to  be  a solution  of  (3.29). 
To  show  wn  6 i/J  we  estimate,  using  (3.6) 


(3.32)  «wnlro  < IlgJI^k  l e-ot-*k  = H g U 

j = 0 n 


K 


1-  e 


-ak 


g II  K 


(using  x/(l-e  ) ->  1 as  x -*■  0)  and  this  yields  (3.30)  after 
a redefinition  of  K. 

To  show  uniqueness  is  trivial,  if  w were  a bounded 
solution  to  the  homogeneous  equation  (3.2)  we  would  have 
for  n 21  j , 


( 3. 33)  w = W W-1  w . 

n n 2 D 

and  wn  = 0 follows  immediately  upon  letting  j ->  -a.. 

It  remains  to  show  that  (An,An)  is  admissible.  In 

fact  if  N is  a common  e-almost  period  for  C and  q 

n 

we  would  have  wn+N~  wn  to  be  the  unique  solution 

to  an  inhomogeneous  equation  with  inhomogeneous  term 


(3.  34) 


k [Cn+N  Cn]wn+N 


*n+N 


n 


and  for  fixed  k > 0 the  almost  periodicity  of  w is  a 

n 

consequence  of  (3.30). 

We  can  remove  the  factor  1/k  in  the  denominator  by  using 
the  expression  (2.11)  for  Cn  , 


(3. 35) 


Cn  C(kAn' ’ ’ * ,kAn+k)  ' 


where  C is  a smooth  function  of  its  arguments.  Expanding 
Cn+N-Cn  usi-n9  (3.35)  we  see  that  if  N is  a common  e-almost 
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period  for  A^  and  f^,  then  N will  be  an  0(e) -almost  period 
of  , with  the  bound  in  the  0(e)  term  independent  of  k. 

Note  that  these  two  paragraphs  follow  only  from  (i/^L^) 
admissibility  of  (3.29)  and  do  not  depend  on  the  stability 
of  (0.1). 

The  convergence  of  wn  to  the  solution  yn  also  follows 
very  easily.  If  y is  the  unique  AP  solution  to  (0.2)  we 
have 


(3.  36)  [ a.y  , .-k  F S.y 

j Tn+j  4 "Vn+D 


= kf  l <*.j 
j J J 


(y(nk+0kj)  - y(nk))  d0 


~ I 3 ■ ty  (nk+k  j)  - y(nk)]]  = k g(k,nk)  . 

j J 


(The  summation  is  from  0 to  Z in  the  above.)  Here  g(k,t)  is 
AP  in  t for  each  k and 


(3.37)  II  g (k  , t ) II  ^ = o(l)  , k 0 ; 

(3.37)  follows  from  the  uniform  continuity  of  AP  functions  on 
the  whole  real  axis. 

then  (3.36)  implies  that  w is  the 

n 

. n 

unique  AP  solution  to  the  same  equation  as  w^  except  for 
an  error, 


Now  if  w = 
n 


n+£-l 

i 


(3.  38) 


f g(k,nk) 
0 
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and  (3.37)  together  with  (3.30)  yield 


(3. 39) 


11  wn'wn"  = o(l)  , k + 0 . 


We  can  obtain  finer  convergence  results  by  imposing  some 
smoothness  conditions  on  A(t)  and  f(t).  In  fact  if  the  scheme 
has  order  of  accuracy  p,  and  A and  f have  p+1  derivatives, 
then  y also  has  p+1  derivatives.  It  is  shown  in  Henrici 
(p.  247)  that 


(3.40) 


l Vn+:fk  l Vn+j"  i k?+1 


Gil  y 


P+1.1 


where  the  constant  G depends  only  on  the  scheme.  (3.40) 
together  with  (3.30)  yield 


(3.41)  llwn-^nBoo  = 0(kP)  * 

We  point  out  that  these  convergence  arguments  are  valid  when- 
ever we  have  admissibility  and  do  net  require  stability. 

Before  leaving  the  stable  case  we  would  like  to  discuss 

problems  relating  to  the  computability  of  the  solution  w . 

n 

Consider  first  errors  in  initial  data.  Suppose  we 

solve  the  exact  difference  equation  (3.29)  but  use  as  initial 

data  Wq  (k)  = wq  OO  + e (k)  where  Wg(k)  is  the  exact  initial 

data  for  the  AP  sequence  solution  wn  and  e (k)  is  bounded  for 

small  k,  i.e.  He(k)ll  < en.  Now  the  solution  w*  , which  we 

w n 

solve  for,  will  be 


(3.42)  w*  = w + W e (k)  , 

n n n 

and  we  have 
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(3.  43) 


n > 0 . 


Iwn"wnl  = 0(e~akn)  / 

Now  in  general  we  cannot  expect  e(k)  + 0 as  k + 0. 

This  would  certainly  occur  if  we  could  obtain  initial  data 
consistent  with  the  initial  data  of  the  solution  y to  the  ODE. 
However  this  solution  is  known  only  by  an  improper  integral 
involving  not  only  the  given  forcing  function  f(s)  but  also 
the  fundamental  solution  Y(t)  to  the  homogeneous  system 
which  we  could  not  expect  to  know  explicitly  unless  A were 
a constant.  Thus  the  most  we  can  assert  is  that,  using  fixed 

■ • t k 

initial  data  (say  Wg(k)  =0),  if  we  integrate  over  a sufficiently 
long  interval  we  will  get  the  solution  wn  up  to  some  exponen- 
tially decaying  error. 

Now  let  us  consider  round-off  error.  If  we  could 
postulate  that  we  solved  a perturbed  equation 

(3. 44)  wn+1  = CRwn  + kfn  + 0(k  1)  , n > 0,  r > 0 , 

where  the  error  is  uniform  in  n then  we  could  relate  this 
* 

solution  to  w by 
n 1 

(3.45)  w**  = w*  + k l W w”1  0(kr)  , n > 0 , 

j = l 3 

and  we  can  estimate  the  error  simply, 
n , °° 

(3.46)  Ilk  l W W~  0(kr)H  < krKk  Y e_0lkr  = 0(kr)  . 

j=l  n 3 r=u 

Let  us  further  examine  the  validity  of  (3.44).  If  we 
include  rounding  errors  in  forming  the  companion  matrix  (or 
operator) , together  with  the  error  of  taking  only  a finite 
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expansion  of  the  term  (ot£I  - g£kAn+£)-1  in  the  implicit 
case  (see  (2.11)),  then  we  could  assume  that  the  linear 
part  of  (3.44)  would  be  replaced  by 


(3.47) 


Cn  = Cn  + 0(k  > • 


If  we  also  assume  an  error  0 (kr+1ll  wnll  ) in  forming  the 
product  cnwn  we  see  we  should  postulate  a system 

( 3. 48)  w**i  = Cw**  + kgn  + 0(kr+1)  + 0(kr+1lwn**ll) 

and  in  order  to  justify  (3.44)  we  must  show  that  the  solution 
to  (3.48)  will  be  bounded  for  n >_  0. 

We  can  write  this  solution  as 

n 


(3.49) 


* * * * p -.1 

wn  - Vo  + k £ wnwj  9j-l 
n 

+ k[  l W W"1(0(kr)  + 0(kr8w**  1)1 
j=l  n 3 3-1 


Now  if  we  use  the  fact  that 

(3.50)  k l e'“kr  = 0(1) 

r=0 

and  define  h = max  II w.  II  we  can  estimate  h bv 
n j=0 , . . . , n 3 n 

(3.51)  hn  = 0(1)  + kr  hn  0(1) 

where  the  0(1)  terms  are  independent  of  n This  implies, 

small  k/  the  boundedness  of  II  w^  II  and  so  justifies  (3.44). 
The  techniques  used  in  Section  2 of  block  diagonalizating 
the  companion  matrix  were  first  developed  by  Engquist 
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although  I was  not  aware  of  his  work  when  my  formulation  was 
developed. 

Engquist's  procedure  may  appear  to  be  more  complicated 

than  the  procedure  presented  here,  but  this  is  because  he 

proves  a more  general  result;  namely  that  the  companion 

matrix  can  be  block  diagonalized  to  within  an  error  of 

Pi 

arbitrary  order  (say  0(k  x)  ) . (He  also  does  not  restrict 
himself  to  strongly  stable  schemes,  but  this  imposes  stability 
requirements  on  several  additional  homogeneous  systems.) 

If  one  were  to  consider  his  procedure  restricted  to  a block 
diagonalization  up  to  0(k  ) then,  while  the  two  formulations 
differ  in  the  lines  of  approach,  they  are  equally  simple. 

Engquist  also  shows  that  the  upper  block  will  agree 
with  the  Taylor  series  expansion  of  (0.1),  up  to  terms  of 
order  p,  if  p < p^^  is  the  accuracy  of  the  scheme.  The 
proof  of  the  first  order  structure  of  the  upper  block  (the 
term  D(z)  in  Section  2)  is,  I believe,  somewhat  simpler 
than  his  proof  restricted  to  the  first  order  term. 

I would  like  to  point  out  that  in  the  case  we  are 
considering  (an  exponentially  stable  ODE  and  a strongly  stable 
scheme)  the  extra  fineness  of  the  block  diagonalization  being 
carried  to  order  p+1  will  only  give  the  advantage  of  the 
homogeneous  system  having  a stability  exponent  -a+0(k^) 
instead  of  -a+0(k).  In  particular  in  Engquist's  study  of 
the  uniform  convergence  (for  t >_  0)  of  the  solutions  of  the 
homogeneous  equation  to  the  solution  to  the  ODE  (see  his 

-36- 


Theorem  3,  p.  24-27  and  note  that  what  he  uses  as  a is 
what  we  call  -a)  one  would  still  obtain  uniform  convergence 
of  0{kp)  (neglecting  round-off  and  errors  in  initial  data) , 
however  the  decay  exponent  would  be  -a+0(k)  instead  of 
-a+0  (kP)  . 
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4 • Admissibility  in  the  Case  of  a 
General  Exponential  Dichotomy 


Here  we  are  concerned  with  admissibility  properties 
in  the  case  that  the  homogeneous  equation  (0.1)  has  a 
general  exponential  dichotomy  as  described  by  (0.6)  with 
P2  ¥ O’  The  results  will  be  valid  in  the  infinite  dimen- 
sional case,  although  in  this  case  one  must  bear  in  mind 
that  the  situation  described  by  equations  (0.6)  is  not 
the  most  general  form  of  an  exponential  dichotomy. 

First  of  all  consider  arbitrary  homogeneous  and 
inhomogeneous  difference  equations 


(4.1) 

(4.2) 


w 


n , i = C (k)w 
n+1  n n 


w 


n+1  ' Cn(k)wn  + kgn 


Here  Cn  is  defined  for  k e (0,kQ]  and  is  AP  in  n for 


fixed  k,  and  is  now  assumed  invertible  so  that  a fundamental 


solution  operator  Wn  , as  described  by  (3.3),  exists  and 


is  invertible  for  all  n. 

Associated  with  (4.2)  we  have  variation  of  constants 
formulas  in  both  the  forward  and  backward  directions, 

n 


(a) 


w = W w 
n n 


n + k y W W71q.  _ 
0 j = i n 3 y3-l 


(4.5) 


n > 0 


(b) 


0 


Wn  = Wnwn  ~ k l w wT1g  . . , 
0 j=n+l  n J 3"1 


n < 0 


as  can  be  easily  verified. 
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Now  we  will  define  an  exponential  dichotomy  for  (4.1) 
as  follows.  Suppose  there  exist  projections  P^k),  i>2(k) 

= I “ such  that  the  estimates 

"ViW'1'  1 I<  e"ak(n"j)  , n > j , 

(4.6) 

||wnP2W-1|i  < K e-ak  ^ D-n)  , j > n , 


hold  with  K and  a independent  of  k.  The  only  difference 
between  this  and  the  ODE  case  is  that  the  projections 
may  depend  on  k. 


If  we  have  an  exponential  dichotomy  then  (Ii£,Ln)  is 
admissible  for  (4.2).  In  fact  the  unique  bounded  solution 

w is 
n 

(4- 7)  wn(k)  = k l <£g._ . 

j=  — 00  J J **■ 

where 


(4.8) 


Gn  = W P.wT1  , 
J n 1 J 


n 1 j 


= -w  P-wT1, 

n 2 j 1 


j ^ n+1  . 


Equation  (4.7)  is  of  course  suggested  by  the  ODE  case  and 
is  easily  verified  using  (4.6).  The  admissibility  bound 
(3.30)  can  be  derived  exactly  as  in  the  stable  case  as  can 
(An,An)  admissibility. 

We  next  observe  that  (l£,l£)  (or  (An,An))  admissibility 
is  preserved  under  perturbations  of  the  linear  term.  In  fact 
if  we  had  a system 
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U-9)  wn+l  = Cn(k)wn  + kSn(k)wn  + k^n 

where,  for  small  k, 


(4.10)  II S (k ) II  < 

n oo  — o 

and  eQ  is  some  number  to  be  determined,  then  the  unique  Ln 
solution  can  be  found  by  defining  w = 0 and  w . to  be 
the  unique  Ln  solution  to 

OO 


(4.11) 


w 


n+1;  i+1 


C w . , . + kS  w . + kq 
n n , l+l  n n,i  yn 


(for  simplicity  the  dependence  of  C and  S on  k has  been 

n n 

suppressed) . If  < 1,  where  is  the  admissibility 

bound  for  (4.2),  the  contracting  mapping  principle  establishes 
admissibility  with  admissibility  bound 


(4.12) 


K. 


K2  (1-£oK1) 


Of  course  we  will  get  AP  solutions  if  D +kS  and  q are  AP . 

n n ^n 

Now  these  two  principles  certainly  settle  the  question 

of  admissibility  when  the  homogeneous  system  has  a general 

exponential  dichotomy.  Referring  to  the  reduction  obtained 

in  Section  3 we  see  that  after  a nonsingular  change  of 

dependent  variable,  v = T(kA  )w  , v satisfies  the 

n n n n 

equation  (see  (3.8)  ff.); 

r I+D  (k)  0 

( 4 • 1 3)  v = I n vo  + k 0 ( 1)  v + kg 

n+i  l o B J n n 

n 

Here  gn  = T(kAn+^)gn>  By  the  boundedness  of  T and  T~^, 


-40- 


admissibility  for  (4.13)  is  equivalent  to  admissibility 
for  (4.2).  But  (4.13)  expresses  the  homogeneous  term  as 
a perturbation  of  the  system 


(4.14) 


X 


n+1 


I+D  0 
n 

0 B 


X 


n 


n 


with  fundamental  solution 


X 


n 


t Y 0 

n . 
n-1 

t o FT 


(4.15) 


X 


-n 


FT  b.  > 

i=l  1 


r Y 0 ^ 

1 -n 

-n+1  . 

o rr  b-1] 

i = 0 


n > 0 , 


n > 0 , 


and  this  has  an  exponential  dichotomy  with  projections 
independent  of  k and  in  fact  (after  identifying  with  E) 


(4.16) 


Px  = © P(0) 


P = P 
2 v 2 


Here  P^  and  P2  are  defined  in  (0.6)  while  P(0)  is  the 
projection  onto  E^  along  E,  (see  (2.19)  ff) . 

We  thus  obtain  admissibility  in  this  case  and  we  point 
out  that  the  convergence  proofs  given  in  Section  3 did  not 
require  stability  and  remain  valid  here. 


5 • Existence  of  an  Exponential  Dichotomy  for 
the  Homogeneous  Linear  Equation 


We  will  prove  here  that  admissibility  for  the  inhomo- 
geneous difference  equation 

(5*1)  w = C w + kg 

n+1  n n yn 

implies  an  exponential  dichotomy  for  the  homogeneous  equa- 
tion 


(5.2) 


w . . = C w 
n+1  n n 


The  proof  is  very  similar  to  a proof  of  the  analogous  property 
for  the  ODE  case  given  in  the  Appendix.  It  is  a strictly 
finite  dimensional  proof  and  it  may  then  be  simpler  to  think 
in  terms  of  matrices  rather  than  linear  operators.  We  can 
then  regard  T(z)  as  a similarity  transformation  which  follows 
an  initial  transformation  putting  U(0)  into  block  diagonal 
form.  The  matrix  L(z)  is  strictly  in  block  diagonal  form 
and  it  is  no  longer  necessary  for  us  to  continually  make  the 
qualification  "identifying  with  E." 


If  we  let  vn  T(kAR)wn  , (5.1)  is  transformed  into  (se< 


(4.13)  ) 


(5.3)  v 


n + 1 


I+kD  0 
n 


v + k 0 ( 1)  v + kg 
n n rn 


n 


(For  simplicity  we  will  write  the  inhomogeneous  term  in  (5.13 


as  g rather  than  g . ) 
n 
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If  we  call  the  leading  matrix  in  (5.3)  Rn (k) , then  the 
unperturbed  linear  system 


(5.4) 


n+1 


R x 
n n 


has  an  exponential  dichotomy  with  projections  P^  and  P2  , 
independent  of  k (see  (4.15)  and  (4.16)) , and  we  must  exhibit 
an  exponential  dichotomy  for  the  homogeneous  version  of  (5.3). 


(5.5)  vn+1  = Rnvn  + k20(l)vn  . 

By  admissibility  (5.5)  can  have  no  solution.  Let 

S (k)  and  S (k)  be  the  subspaces  of  initial  data  which  give 
rise  to  solutions  bounded  for  positive  and  negative  n respec- 
tively. We  have  S+  n S_  = {0}  (suppresssing  the  k dependence) . 
If  we  let  S be  any  complementary  space  we  have 


(5.6)  e£  = s+  e s_  e s 

with  associated  projections 


(5.7)  I = P1  + P2  + P3 

Our  first  task  is  to  show  S = {0}  i.e.  P^  - 0. 

^ti  n - 

Define  to  be  the  submanifold  of  Lot  of  sequences 

with  only  finitely  many  nonzero  components.  Let  gn  e L^ 

and  let  vn  be  the  corresponding  unique  L^  solution  to 

(5.3).  Then  using  (4.5a)  we  see  that  for  large  positive  n 

we  have 

00 

(5.8)  vn  = Vv0  + k ^ Vj  Sj-l1 
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is  bounded  the  vector  in  the  brackets  lies 


and  since  v 


n 


in  S+  i.e, 


(a) 


,v  = - k l P~  V_1  g.  . 
° -j  = i 2 3 


(5.9) 


(b) 


P3V0  * - k J,  P3  Vj1  9-j-l 
3 = 1 J J 


Similarly  looking  at  vn  for  n and  using  (4.5b) 

obtain 


we 


(a) 


lv0  = k . I P1  Vj1  9-j-l 

j = — 00  J J 


(5.10) 


(b) 


P3V0 


= k l P , vT 1 g . 

. 3 ■)  ^i-l 

J = — oo  -J  j 


Now  if  we  simply  set 


■ V1P3Z 


(5.11) 


9^  — 0 / i ^ 0 , 


where  z is  an  arbitrary  vector  in  we  see  immediately 
that  P^  = 0.  We  can  also  see  that  for  any  given  g 6 


. n 


the  unique  Lro  solution  vn  is  given  by 


(5.12) 


v 


• = k l Gn  g.  . 

n jf-00  3 y3"l 


where  is  given  by  (see  (4.8)) 


(5.13) 


Gi  “ Vn  P1  Vi1  ' n - 3 


."I 


“ _Vn  P2  V j ' 5 1 n+1  • 
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It  is  of  course  which  we  wish  to  prove  is  exponentially 

damped. 

Let  Gj  be  the  Green's  function  associated  with  (5.4). 
We  know 

(5.14)  IlGjll  < k e"akln_:i  I 

Now  define  Un  = G1}  - G1}  . If  g.  e and  v the  corres- 

J J J J 00  n 

ponding  solution  to  (5.3),  while  x is  the  Ln  solution 

n oo 

to  the  inhomogeneous  version  of  (5.4), 


(5.15) 


Then  u 

n 


(5.16) 


v 

n 


x , , = R x + kg 
n+1  n n yn 

- x^  is  the  unique  bounded  solution  to 

u . -i  = R u + k2  0(1)  v . 
n+l  n n n 


We  thus  have 


(5.17)  u = k l Hn  g . . = k2 

n j=-oo  3 yD-l 

If  we  fix  i and  let  g.  . = z, 

^i-l 

z is  an  arbitrary  vector  in  , 
and  the  admissibility  bound  for 


l Gn  0(1)  v . 

jz=  — 00  J J 

gj  = 0 for  j ^ i-1  where 
we  then  obtain,  using  (5.14) 
(5.3)  , 


(5.18)  « Hj  II  < K1  , for  all  n,  i 

where  is  some  constant  independent  of  k.  (The  restric- 

tion "for  small  k"  is  always  understood.)  The  use  of  (5.14) 

then  implies  the  boundedness  of  Gn  i.e. 

3 

(5.19)  II G1? II  < K , 

D ^ 

where  K ^ is  independent  of  k. 
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We  will  now  show  that 


(5.20)  II  P 1 ( k ) - = 0 (k)  . 

Observe  that  P^  is  independent  of  k.  To  obtain  (5.20)  we 
note  that  P^-P^  = we  define  = z'  9j  = ® 

(j  ¥ i)  where  z is  an  arbitrary  vector  we  then  have 
from  (5.12) 

(5.21)  v.  , = k . 1 g.  . = k G?  1 z . 

j-1  i 3i-l  i 

Substituting  into  (5.17)  we  obtain 


z 

Now  using  (5.19)  and  (5.14)  we  see  the  sum  in  brackets  in 

(5.22)  is  0 ( 1/k ) whence  we  conclude  that 

(5.23)  iHjl  = 0 (k) 


(5.22)  u = kH*?z  = k2 
n i 


3=-° 


5"  0(1)  G^_1 


and  setting  n = i = 0 we  obtain  (5.20). 

Note  that  this  proves  that  the  ranks  of  the  stable  and 
unstable  manifolds  are  unchanged  for  small  k. 

We  can  now  show  quite  easily  that  (5.19)  can  be  replaced 
by  an  exponential  decay  factor. 

For  a certain  small  positive  e consider  new  systems 

• • ™ 
with  the  linear  part  of  (5.3)  and  (5.4)  multiplied  by  e i.e. 


(5.24)  v^+1  = e ck[Rn  + k2  0(1)  ] v^ 

(5.25)  x,.=e  Rx 

n+1  n n 
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Now  by  Section  4 (see  (4.9)  ff.)  the  inhomogeneous  versions 
of  (5.24)  and  (5.25)  will  both  have  (Ln,Ln)  admissible,  if 
e is  small  enough.  But  the  new  fundamental  solutions  are 

(5.26)  v£  = e_enk  V 

n n 

(5.27)  XE  = e'£nk  X 

n n 

and  clearly  (5.25)  will  have  an  exponential  dichotomy  with 
projections  P1  and  P2  if  e is  small  enough. 

Now  if  P-^(k)  and  P2  (k)  are  the  corresponding  stable 
and  unstable  projections  for  (5.24)  while  G1?  # is  the 
Green's  function  for  (5.24),  the  same  analysis  that  led  to 
(5.19)  will  yield  for  some  constant  K ^ , 

(5.28)  II  Gn  II  < K 0 

D ; e - 3 

In  particular  if  j n+1,  (5.28)  reduces  to 

(5.29)  II V P£  vT1!  < K,  e“ek(^_n) 

n 2 j - 3 

and  one  part  of  the  definition  of  an  exponential  dichotomy 

p 

for  (5.5)  will  have  been  shown  if  we  can  show  P^  = P^  , 

E 

P2  = P2  . But  this  is  a simple  consequence  of  (5.20).  In 

EE 

fact  if  S+  , S_  are  the  new  stable  and  unstable  manifolds 
for  (5.24)  then  (5.26)  yields 

(a)  S,  c s£ 

(5.30)  + 

(b)  sf  C s_ 

E 

But  (5.20)  shows  that  the  ranks  of  P^  and  P^  are  equal 
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(being  equal  to  the  rank  of  P-^  and  similarly  for  and  P2 . 

This  shows  that  the  inclusions  in  (5.30)  are  equalities  and 
that 

(a)  P^  = P, 

(5.31)  E 

(b)  P^  = P2 

Thus  one  part  of  the  requirements  of  an  exponential  dichotomy 
for  (5.5)  has  been  shown  and  the  other  part  will  follow  on 
replacing  -e  by  +e. 

Lastly  we  point  out  that  an  exponential  dichotomy  for  (5.5) 
implies  one  for  (5.1)  since  in  terms  of  fundamental  solutions 
we  have 

(5.32)  Wn  = T-1(kAn)Vn  T(kA0) 

whence  the  projections  P^(k),  P2  (k)  for  (5.2)  are  related  to 
p^  and  P2  by  a similarity  transformation 

(a)  P^(k)  = T"1(kA0)P1(k)T(kA0) 

(5.33) 

(b)  P2(k)  = T_1(kA0)P2(k)T(kA0)  . 
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6.  Miscellaneous  Results 


This  section  is  concerned  with  certain  generalizations 
of  the  theory  developed  in  the  preceding  section.  We  first 
remove  certain  restrictions  which  had  been  imposed  previously 
in  order  to  make  the  exposition  clearer.  In  Section  6A  we 
will  remove  the  restriction  that  A be  while  in  Section  6B 
we  remove  the  restrictions  (2.6)  and  (2.7)  on  the  roots  of 
the  polynomial  p(x)  = 0 which  lie  inside  the  unit  circle. 

We  then  consider  certain  trivial  extensions  of  the 
theory.  In  Section  6C  we  remark  on  the  general  (non-AP) 
case.  In  Section  6D  we  consider  the  convergence  of  the 
mean  value  of  our  sequence  solution  to  the  mean  value  of 
the  solution  of  the  ODE  and  finally  in  Section  6E  we 
consider  the  weakly  nonlinear  case. 


; 


6A*  — moval  of  the  Differentiability  Condition  on  a 


First  consider  the  case  when  the  homogeneous  equation 


(6.1) 


y = Ay 


is  exponentially  stable;  i.e. 


(6.2) 


Y ( t)  Y-1  (s)  11  < K e_a(t_s)  , t 


> s. 


where  Y(t)  is  the  fundamental  solution  to  (6.1).  We  assume 
the  A is  AP  but  is  not  C1.  Now  Theorem  2 does  not  require 


a perturbation  0(k  ) but  is  certainly  valid  for  a perturbation 


k o(l)  (k  -*■  0).  This  leads  us  to  expect  that  the  requirement 
.1 


that  A be  C can  be  replaced  by  the  uniform  continuity  of  AP 
functions  on  the  whole  real  axis.  This  is  indeed  the  case 
as  one  can  verify  with  little  difficulty. 

In  this  subsection,  we  will  merely  trace  through  the 


proofs  in  Section  3 and  indicate  what  changes  must  be  made 
.1 


if  A is  not  C . 

Discretizing  (6.1)  we  obtain  the  homogeneous  difference 
equation 


(6.3) 


W , = C w 
n+1  n n 


where  Cn  is  given  in  (2.11).  Using  the  uniform  continuity  of 


A we  see  that  (3.7)  can  be  replaced  by 


(6.4) 


wn+l  ~ U(k  An)wn  + k o(1)wn 


where  the  term  o(l)  is  uniform  in  n (this  will  not  be  stated 
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explicitly  in  the  future) . If  we  also  use  the  uniform 
continuity  of  A to  obtain 


(6.5) 


T(k  An+lJ  = T(k  V + k o(1) 


we  see  that  Theorem  2 will  yield  exponential  stability  for 
(6.3)  provided  we  can  show  exponential  stability  for  the 
system 


(6.6) 


v 


n+1 


I+kA  0 
n 


B 


v 


n 


n 


The  lower  block  causes  no  difficulty  and  we  are  left  with 
the  system 


(6.7) 


x11=  [I  + k A ] x . 
n+1  n n 


Finally  we  can  show  that  (6.7)  is  exponentially  stable  by 

using  Theorem  2 to  compare  the  fundamental  solution  X with 

n 

Yn  (=  Y (nk) ) . This  reauires  some  modification,  as  the  proof 
leading  to  (3.16)  used  the  differentiability  of  A (see 
(3.13)  ff ) . 


Following  the  notation  of  Section  3,  we  define 

-1 


H ( t , s ) = Y ( t ) Y (s) . We  can  then  write  (compare  with  (3.14)) 

1 

(6.8)  H ( (n+1) k,  nk)  = I + kH(nk,nk)  + k j d0 [H (nk+0k ,nk) 


- H (nk ,nk) ] 


= I + kAn+  k d6  [A (nk+8k) -AnJ  H(nk+6k,nk) 


+ k d6  An [H (nk+0k,nk)  - H(nk,nk)] 
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Finally  we  see  that  if  we  use  (6.2)  together  with  the 
uniform  continuity  of  A on  the  whole  real  axis  and  the 
equation 

(6.9)  H ( t , s ) = A ( t)  H ( t) 
then  (6.8)  can  be  rewritten  as 

(6.10)  Yn+1  = [I  + kAn  + k o(l)]  Yn 

and  Theorem  2 is  immediately  applicable. 

We  thus  have  exponential  stability  for  (6.3)  and  hence 
(An,An)  admissibility  for  the  inhomogeneous  version  of  (6.3). 

Finally  we  point  out  that  the  results  of  Section  4, 
regarding  admissibility  in  the  case  that  (6.1)  has  a general 
exponential  dichotomy  is  equally  valid  if  A is  not  C1,  as 
the  contracting  mapping  principle  (see  (4.9)  ff)  would 
certainly  be  applicable  if  the  perturbation  is  k o(l).  The 
results  in  Section  5,  however,  use  crucially  the  differenti- 
ability of  A (in  the  argument  involving  the  deduction  of 
(5.18)  from  (5.17))  and  it  has  not  been  possible  to  extend 
this  result  when  A is  not  C^. 
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1 


6B.  Removal  of  Restrictions  on  the  Roots  of  p(x)  = 0 

In  Section  2 we  prescribed  two  conditions  on  the  roots 
of  p(x)  = 0 which  lie  inside  the  unit  circle;  namely 

(6.11)  xu  distinct  , u = 2, ...,£, 

(6.12)  xu  ¥ ° , u = 2,  . . , 1. 

(Recall  that  x^  = 1 and  | xu | < 6 <1  for  u > 1) . We  will 

now  remove  these  restrictions. 

It  should  be  clear  by  now  that  in  the  exponentially 

stable  case  (6.12)  is  unnecessary.  In  fact  the  fundamental 

solution  to  the  homogeneous  equation  is  used  only  for  n ■>  j , 

and  the  equation  (3.5)  shows  that  the  invertibili ty  of  the 

linear  term  is  not  required  for  n ^ j.  Of  course  the 

notation  W W.^  is  no  longer  accurate,  but  except  for  this 
n ] 

detail  the  results  of  Sections  2 and  3 are  valid  without 
assuming  (6.12) . 

We  deal  next  with  the  restriction  (6.11).  (Note  that 
the  root  x^  = 1 is  always  simple.)  The  removal  of  (6.11) 
in  the  finite  dimensional  case  is  trivial.  In  fact  suppose 
the  root  xr  has  multiplicity  q > 1;  i.e. 

(6.13)  P1(xr)  = ...  = pq"1(xr)  = 0 
Pq(xr)  f 0 . 

It  is  shown  in  Henrici  (p.  214)  that  the  effect  of  (6.13) 
is  that  the  eigenvalue  x^.  of  U(0)  (see  (2.12))  will  now  have 
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nontrivial  Jordan  blocks.  Thus  U(0)  can  still  be  put  in 
block  diagonal  form  (see  (2.17)) 

f 1 0 

(6.14)  U ( 0 ) = 

1 0 B (0) 

where  B(0)  is  no  longer  diagonalizable , but  has  all  of  its 
eigenvalues  bounded  by  6 in  the  norm.  Now  since  B(0)  can 
be  put  in  Jordan  normal  form  with  e instead  of  1 on  the 
superdiagonal,  for  any  e > 0,  we  can  certainly  choose 
a matrix  norm  such  that 

(6.15)  II  B ( 0 ) U < 0 


Since  all  norms  are  equivalent  on  a finite  dimensional  space, 
the  proof  of  Theorem  1 can  now  proceed  exactly  as  in  Section  2. 

The  case  when  E is  infinite  dimensional  can  be  handled 
in  exactly  the  same  manner?  however  since  we  can  no  longer 
appeal  to  the  theorems  of  linear  algebra,  we  will  have  to 
carry  out  the  proof  in  more  detail. 

It  is  shown  in  Henrici  (p.  214) , that  in  the  scalar  case 
a basis  for  the  generalized  eigenspace  Er  corresponding  to 
the  eigenvalue  x is  {e  . } where 


(6.16) 


f X*"1  1 
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( j-l 


fT  U-r)x4'  ‘I'11  ] 

r=l 


ff  (n-r)x(n_(j_1)) 
r=l 


fr  (j  - r)x 


r=l 


j > 1 


One  can  verify  easily  that  in  fact 


(a)  U(O)  ir>1  = xrlrjl 


(6.17) 


(b)  U (0)  e . = x e . + e . . , 
r,D  r r,D  r,j-l 


3 > 1, 


Also  since  the  vectors  e . form  a new  basis  the  l x l 

r > 3 


matrix  with  columns  e . is  nonsingular.  Now  in  order 

r ' 3 


to  have  e instead  of  1 on  the  superdiagonal,  we  replace  e 


r,D 


by  e . where 


(6.18) 


j-l  - 

e . = eJ  e 

r,D 


(see  Bellman  [9],  p.  198).  Here  e > 0 is  to  be  specified. 


The  matrix  with  columns  e . is  still  nonsingular  but 

r,D  y 


(6.17b)  will  now  be  replaced  by 


(6.19) 


U(0)e  .=x  e . + e e ... 

r,D  r r,n  r,]-l 
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Now  the  extension  of  this  to  infinite  dimensional  space 

is  very  simple  and  follows  closely  the  procedure  of  Section  2. 

For  any  vector  y £ E define  the  vector  e . (y)  £ E0  by 

r / ] x, 


(5.20) 


e -(y)  = 1 

r • J 


' FT  U-r)xl-jy  ] 
r=l 


] f ( j + l-r)  x y 
r=l 


Equation  (6.19)  clearly  generalizes  to 


(6.21) 


U(0)  e .(y)  = x e . (y)  + e e . . (y)  , i > 1 

r,]  r r,j  1 r,;j-l  1 ' J 


U(0)  era(y)  = xr  era(y) 


We  now  define  the  space  E . = {w  £ E0  such  that 

r i j x, 

w = e . (y)  for  some  y £ E}.  Clearly  the  space  E . is 
r > J r I 3 

closed  for  all  r and  j.  Furthermore  we  have 


(6.22) 


E.  = < 

a u . r,  i 

r,D  J 


l E ■ 

■ a rO 


This  follows  by  the  same  argument  as  that  given  in  Section  2 

( Xl 

(see  (2.14)  ff)  . In  fact  if  w = ■ 


and  we  had  an  expansion 


j 1 

^ H 


is  a vector  in  E 


£ 


(6.23) 


r f1 

- H 


= J e . (y  . ) 

L>  V . -l  Y-  -| 


r / 3 


then  the  {y  .}  would  be  related  to  the  {x.}  by  multiplication 
r > D 1 
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kMMi. 


■ ‘■■.on -- ...  ■. 


by  a nonsingular  £ x £ matrix  and  this  establishes  (6.22). 

It  also  follows  from  this  that  if  we  define  a new  norm  on 

E.  by  II  wll  u = max  II  y .11  then  II  wll  u is  equivalent  to  II  wll  ^ 

* r-  -■  r'  3 E 


( max  II  x . II  ) . 
i 1 E 


Now  if  is  defined  as  E^  ^ (the  root  x^  is  simple) 


while  E is  defined  as 


I 

r>l,  j 


» then,  exactly  as 


• J 

in  Section  2,  we  see  that  U(0)  is  in  block  diagonal  form 
with  respect  to  the  decomposition  E = E^  © E1,  and  we  can 
write  (see  (2.17)) 


f D (0)  0 

(6.24)  U (0 ) = 

^ 0 B (0) 

where  D(0):  E^  -*■  E^  and  is  the  identity,  while  B(0):  E^"  -»■  E 
and  for  e sufficiently  small  we  have  from  (6.19) 

(6.25)  II  B ( 0 ) II  U < 0 


The  proof  of  Theorem  1 can  now  be  carried  out  exactly 

as  in  Section  2 if  we  make  the  final  observation  that  the 

spectrum  of  U(0)  is  exactly  {x^}.  This  should  be  obvious 

and  follows  from  the  fact  that  on  each  "generalized  eigen- 

space"  E = © l E . , U(0)  - AI  (for  A f {x  })  acts  as 
r j r , j u 

the  matrix 


xr-  X z 


• • 


(6.26) 


0 x -A 
r 


and  the  inverse  is  obtained  merely  by  inverting  this  matrix, 


6C.  Extension  to  the  L Case 


Here  we  would  like  to  make  the  simple  observation  that 
the  difference  equation  theory,  like  the  ODE  theory,  is 
essentially  an  theory  which  produces  AP  solutions  when 
the  coefficients  are  AP. 

In  fact,  if  we  now  agree  that  derivatives  are  to  be 
understood  as  instead  of  AP,  the  theory  developed  in 
Sections  3,  4 and  5 is  entirely  valid  in  the  Lot  case,  with 
only  one  unimportant  exception.  The  single  exception  is 
the  proof  given  in  (3.36)  ff  of  the  convergence  of  the  Ln 

yn+o_i) 

solution  w to  the  vector  w = . in  the  case 

XI  XI  • 

yn 

that  the  coefficients  are  not  smooth.  In  fact  if  we 
refer  to  (3.36)  and  (3.37)  we  see  that  (3.37)  need  not 
be  valid  because  the  function  y which  is  (we  take  A and 
f continuous)  need  not  be  uniformly  continuous  on  the 
whole  real  axis.  In  this  case  we  have  not  been  able  to  show 
that  wn  -*■  wfi  uniformly  for  all  n but  only  uniformly  for  nk 
lying  in  compact  intervals. 


To  see  this  let  us  write  w (k)  , w (k)  to  indicate 

n n 


explicitly  the  dependence  on  k.  Now  if  the  initial  data 


for  w„ (k)  is  consistent,  that  is 
n 


i y- 


( 6.27) 


wo(k)  -w 


l *1 J 


for  some  vector  y^  e E,  then  it  is  well  known  (see 


Henrici,  p.  244)  that  wn(k)  roust  converge  uniformly  on 
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compact  subsets,  to  some  solution  to  the  equation 

(6.28)  y = A ( t)  y + f . 

But  since  every  solution  to  (6.28)  but  one  is  unbounded,  this 
limiting  solution  can  only  be  the  unique  Lro  solution  y(t).  Thus 
we  must  show  (6.27)  and  the  proof  of  this  will  in  fact  show 
that  the  vector  y^  in  (6.27)  is  y(0)  . 

We  assume  at  least  that  A is  (with  bounded  derivative) 
so  that  the  results  of  Section  5 are  applicable.  We  can  then 
write 

00 

(6.29)  wQ  - w = k l G?  g(k,(j"l)k) 

j — oo  J 

where  g is  given  in  (3.37),  and  G^  satisfies  the  estimate 

-ot  k | n-  j | 

(6.30)  H G .11  < K.  e 

3 - 1 

Observe  that  g(k,t)  0 uniformly  for  t in  compact  intervals.  Let 
T be  unspecified  for  the  moment  and  rewrite  the  sum  in  (6.29)  as 

(6.31)  wn-w  = k l G°g(k, ( j-l)k)+k  l G^g (k,  ( j-1) k)  . 

| jk | <T  J | jk | >T  ] 

The  first  term  0 because  g(k,t)  ->  0 uniformly  on  the  compact 

-a^T 

interval  [ -T , T ] . The  second  term  is  0(e  ) as  we  can  easily 

see  from  (6.30).  Thus  by  first  choosing  T sufficiently  large 
to  make  the  second  term  small,  and  then  choosing  k small  enough 
to  make  the  first  term  small  we  see  that  IIwq-WqII  •>  0 (k  -*■  0)  and 
this  establishes  the  convergence  of  wn  to  wn  uniformly  on 
compact  intervals. 
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6D . Mean  Value  Properties 


Let  an  e A . It  is  shown  in  Corduneanu  (p.  48)  that 
the  limit 


(6. 32) 


, n+N-1 

lim  n I a.  = m(a  ) 
j=n  J 


exists  and  that  this  limit  is  uniform  and  independent  of  n, 
which  we  henceforth  set  equal  to  zero.  This  is  of  course 
the  exact  analogy  to  the  mean  value  of  an  AP  function 

T 

(6.33) 


lim  - f(t)  dt  = m ( f ) . 

T+co  j 
0 


Now  if  wn  is  our  AP  sequence  solution,  we  would  like  to 

study  the  behavior  of  m(w  ) as  k -*■  0 . Define  w as 

n n 


(6.34) 


w = 
n 


' Yn+£-l  ' 


■n 


where  y is  the  AP  solution  to  the  ODE 


(6.35) 


y = Ay  + f 


Clearly 

(6.36) 


m(wn)  = 


f m<yn)  ) 


^ m (y  ) 


Now  since  (see  (3.39)) 


(6.37) 


lw  -w  I!  = o(l) 
n n “ 1 

k+1 
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we  must  compare  m(yn)  with  m(y).  We  first  claim  for  any 
AP  function  g we  will  have 


(6.38) 


m(9n>  -kTo->  ">(9) 


where  gn  - g(nk).  This  is  obvious  because  (6.38)  holds  if 
g is  a trigonometric  polynomial  and  for  arbitrary  g,  we 
simply  approximate  by  a sequence  of  trigonometric 
polynomials . 

Now  (6.38),  together  with  (6.37)  and  (6.34),  certainly 
yields 


(6.39) 


m(w  ) 


n k+0 


f m(y)  ) 


m(y) 


We  now  concern  ourselves  with  the  rates  of  convergence, 
assuming  smooth  coefficients.  If  A and  f are  CP+1  (so  that 
y is)  we  can  replace  (6.37)  by  (see  (3.41)) 

(6.40)  II  w -w  II  = 0(kP) 

n n °°  ' ' 


where  p is  the  order  of  accuracy  of  the  scheme.  It 
therefore  follows  that 


(6.41) 


m(w  ) 
n 


' m(yn) 
. m(yn> 


+ 0 (kp) 


and  it  is  only  necessary  to  study  the  convergence  of  m(y  ) 
to  m(y)  . In  fact  for  any  AP  function  y cr+1  we  have 
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(6. 42) 


m(yn)  = m(y)  + 0(kr+1) 


Equation  (6.42)  is  an  almost  periodic  analogue  of  a 
theorem  of  Isaacson  and  Keller  [10]  for  periodic  functions 
(p.  340)  and  in  fact  the  proof  is  very  similar,  and  is  a 


simple  consequence  of  Taylor's  theorem.  Letting  t = nk 

n 


and  s e [0,1]  we  can  write 


(6.43)  y(tn+sk)  = yn  + (sk)yj  +...+  — yj  + 0(kr+1) 


where  the  remainder  is  uniform  in  t since  yr+1  is  AP. 
Integrating  (6.43)  from  s = 0 to  s = 1 we  obtain 

fcn+l  1 

(6.44)  | y ( t J dt  = k | y (tn  + sk)  ds 

fcn  0 

r 


= (~r+T)T  ^n  + 0(kr+1)]. 


If  we  sum  (6.44)  from  n = 0 to  N , divide  by  Nk  and  let  N + °° 
we  obtain 


(6.45)  m(y)  = m(yn)  + | m(yj)  + 


TrTTTT  m(^)  + °(kr+1) 


Now  if  r - 0,  i.e.  y is  the  terms  involving  m(y1) 


through  m(yn)  are  missing  and  we  obtain  (6.42).  For  r > 0 
we  assume  (6.42)  holds  for  j = 0,  ...,  r-1.  We  can  then 


1 

write,  as  y is  C , 


(6.46) 


0 = in  (y X)  = m(yj)  + 0(kr) 
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.... .. 


r since 


and  similar  expressions  for  mfy1)  for  i - 2,..., 

e c(r-i)+l  (we  have  used  the  obvious  fact  that  the  mean 
value  of  a derivative  is  zero),  to  obtain  (6.42).  If  we  now 
apply  this  to  (6.41)  we  obtain 

f m(y) 


(6.47)  m(wn)  = 


+ 0(kP)  . 


I m(y)  J 
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6E.  Weakly  Nonlinear  Equations 


Consider  the  equations 

(6.48)  y = Ay 

(6.49)  y = Ay  + f(t) 

(6  • 5°)  y = Ay  + h (t  ,y) 

If  any  pair  of  function  spaces  (B..D)  is  admissible  for 

(6.49)  the  contracting  mapping  principle  provides  a technique 
to  obtain  D-solutions  to  (6.50)  under  certain  conditions  on 
the  nonlinear  term  h(t,y),  the  most  important  being  that 
h has  a small  Lipschitz  constant  (see  Hartman  [11], 

Chapter  12)  . 

The  case  of  (A, A)  admissibility  is  particularly  simple. 
We  require  that  h be  almost  periodic  in  t,  uniformly  for  y 
in  compact  subsets  (see  Hale  [12]  pp.  113  f f)  . This  simply 
insures  that  h(t,x(t))  will  be  AP  for  any  AP  function  x. 

We  also  require  that 

(6.51)  lh(t,x1)  - h(t,x2)l!oo  £ e1llx1  - x2 II 

for  any  x^,x2  E e.  The  restriction  on  is  simply 

(6.52)  < 1 

where  K1  is  the  admissibility  bound  for  (6.49)  (see  (0.3b)). 
The  unique  AP  solution  to  (6.50)  will  simply  be  the  limit 
of  the  iterates 
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(6.53) 


0 


0 

y = 

i+1 

y 

*i+l 

y 


is  the  unique  AP  solution  to 
= Ay 1+1  + . 


This  is  of  course  an  immediate  consequence  of  the  contracting 
mapping  principle. 

The  same  proof  will  yield  AP  sequence  solutions  to  the 
discretized  version  of  (6.50). 


(6’54)  wn+l  = Cnwn  + k 9(nk'k'wn+1'wn) 

Here  the  nonlinear  term  is 

f (a£I-k3£A(t+£k)“1t3£h(t,x£_1) 

l-l 

+ l 3.h(t,x..)] 

(6.55)  g(t,k,w,w)  = I ^-0 


where  w - 


f *1-1 

r 

*£-l 

and 

w = 

■ 

■ 

x„ 

X- 

o 

o 

Observe  that  g is, 


for  fixed  k,  AP  in  t uniformly  for  w and  w in  compact 
subsets,  and  that  g will  have  an  O(e^)  Lipschitz  constant 
which  we  call  e,  i.e. 


(6.56)  II  g(t,k,w2,w2)  - g ( t ,k  ,w1,w1)  II 

e max  [Ilw2-w^ll  , Iw2-w^ll]  . 

Now  if  K is  the  admissibility  bound  for  the  linear  inhomogeneous 
version  of  (6.54)  (see  (3.30))  and  p = eK  < 1 then  (6.54) 


-65- 


will  have  a unique  AP  solution  wn  , which  in  fact  will  be 

the  limit  of  the  iterates  w . defined  by 

n ; 1 J 


(6<57)  wn ; 0 = 0 


n.^+^  is  the  unique  AP  solution  to 
n+l;i+l  ' Cn"n;i+1  + kS (nk 'k 'wn+l; i 'wn; i > 


This  is  an  immediate  consequence  of  the  contracting  mapping 
principle  and  we  note  for  future  reference  that 


(6.58) 


w -w  II  = 0(pm)  . 

n n;m  ® 


Now  the  first  thing  to  consider  is  the  convergence  of 


w to  w defined  as 
n n 


f y 


n+£-l 


yA 


where  y is  the  solution  to 


(6.50)  (assuming  of  course  that  (6.52)  holds).  In  fact 
as  was  seen  in  Section  3 wn  satisfies  (6.54)  up  to  an  error 

pH*  1 

k o(l)  (0(k^  ) if  everything  is  smooth).  Hence  the 

k-+l 

difference  z = w - w is  the  unique  AP  solution  to 
n n n ^ 


(6.59)  z , = C z + k[g(nk,k,w  , , + z ,w  +z  ) 

n+1  n n r n+1  n+1  n n 

- g(nk ,k,wn+1,wn) ] + k o(l)  (0(kp+1)) 

and  we  have 

(6.60)  II  z II  < pll  z II  + o(l)  (0(kp)) 

and  since  M < 1 this  settles  the  question  of  convergence. 

We  now  examine  the  question  of  the  computability  of 
the  solution  in  the  exponentially  stable  case;  that  is  when 
we  have  constants  K and  a such  that 


(6.61) 


n > D, 


II W wT  1ll  < K e“ak  (n“  j ) 
n D - 

where  is  the  fundamental  solution  to  the  homogeneous 

version  of  (6.54). 

To  do  this  we  use  (6.58)  to  reduce  the  problem  to  a 
study  of  the  effect  of  round-off  errors  and  errors  in  initial 
data  on  the  solution  w for  a fixed  m.  This  will  be  a 

n / in 

simple  consequence  of  the  stability  results  of  Section  3 
(see  (3.42)  ff)  for  the  inhomogeneous  equation. 

In  fact  the  results  of  Section  3 show  that  at  the 
first  iteration  we  would  solve  (for  n >_  0)  for  a solution 

•k 

w , which  is  related  to  w , by 
n ; 1 n ; 1 2 

(6.62)  w*;1  = wn;1  + 0(kr)  + 0(e"akn)  . 


We  now  consider  the 
iterate.  Assuming, 
exact  initial  data, 
the  equation 


effect  of  this  error  on  the  second 

for  the  moment,  exact  calculations  and 

* 

we  obtain  a solution  w „ which  solves 

n } £• 


(6.63)  w 
If  we  let 


* 


n+1 ; 2 


= Cnwn ; 2 + k 9 (nk' k 'wn+l ; l'wn; 1>  ’ 

★ 

= Wn.^  “ wn.j_  f°r  i = 1^2  then  we  obtain 


(6.64)  z , = C z + k[g(nk,k,w  ,,  ,+  z n ,w  ,+z  ,) 

n + l?2  n n;2  ^ n+l;l  n+l;l  n;l  n;l 


- gW,k,,n+1.r,n.1i] 

The  last  term  can  be  written  (using  (6.62))  as 
k 0(kr)  + k 0 (e  akn)  and  we  have 
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(6.65)  z 


n;  2 


= k l WwT1  0(kr)  + k 
j=l  3 

= 0(kr)  + 0 (e~akn  kn)  = 


I W W"1  0(e‘akj) 
j = l n 3 

r -a.kn 

0 (k  ) + 0(e  1 ) 


where  can  be  arbitrarily  close  to  a.  We  have  used 

(6.62),  (6.61)  and  the  trivial  fact  that  t = 0(e<5t)  for 

any  positive  6 . 

It  is  now  obvious  that  (6.65)  is  preserved  if  we 

include  in  the  equation  for  w*;2  (see  (6.63))  the  effect 

of  roundoff  and  errors  in  the  initial  data  of  w* 

n;  2 

(compared  with  the  initial  data  for  w ) and  it  is  also 

Xl  / £ 

obvious  that  this  process  can  be  repeated  for  m itera- 
tions, for  any  fixed  finite  m,  and  this  establishes  the 
computability  of  the  solution  w . 


■a.-.—  
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4 





Part  II.  Extensions  to  Partial  Differential  Equations 


We  will  now  use  the  techniques  of  Massera  and  Schaffer, 
in  the  simple  exponentially  stable  case,  to  obtain  AP  solu- 
tions to  a certain  class  of  inhomogeneous  partial  differential 
equations  (PDE's)  to  which  one  can  apply  ODE  formalism  by 
use  of  the  Hi lie  Yoshida  Theorem.  The  homogeneous  equation 
has  been  extensively  studied  by  Krein  [13],  to  which  we  will 
refer  often.  A more  succinct  study  of  the  Hille-Yoshida 
Theorem  can  be  found  in  the  appendix  of  Lax  and  Phillips  [14]. 

Throughout  Part  II  E will  denote  an  infinite  dimensional 
Banach  space. 
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7.  Existence  Theorems 


Consider  the  linear  equation 

(7.1)  y = By 

where  B is  an  unbounded  operator.  The  study  of  the  solutions 
to  (7.1)  has  been  answered  by  the  Hille-Yoshida  Theorem  in 
the  case  that  the  resolvent  R^(B)  (=  (Al-B)  1)  satisfies 

the  inequality 

(7.2)  " ra  (B)  II  < j , A > 0 . 

A more  general  formulation  is  given  in  Krein  (in  particular 
Section  2,  Chapter  1).  Krein  shows  that  if: 


(7.3) 

(a) 

B is  closed 

(b) 

B is  densely  defined 

(c) 

■ r?(B)I  < , R.P.  A > 

(Re  A-  w)n 

to  , 

then 

there 

exists  a semigroup  zt  > t >_  0 , such 

that: 

(7.4) 

(a) 

is  strongly  continuous 

(b) 

tsi 

o 

II 

M 

(c) 

II  Ztll  < Me wt 

(d) 

lim  [ZA-I]e/A  converges  iff  e e Dn 
A-*-0  A B 
and  in  that  case  it  converges  to  Be 

(e) 

if  e £ then  for  t tg  the  unique 

solution  to  (7.1)  such  that  y(tg)  = 

e 

is  simply  y(t)  = Z e (uniqueness 

t-t0 

is 
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shown  in  Theorem  2.7,  page  47) 

(f)  : Dfi  ->  Dfi  and  commutes  with  B for  t > 0 . 

For  a well  posed  problem  we  obviously  require  w < 0 in  (7.4c) 
To  insure  exponential  stability  we  are  going  to  proceed  as 
if  (o  = 0 and  introduce  our  own  damping  term  which  will  be 
allowed  to  depend  on  t.  The  reader  should  note  that  if 
w < 0 this  term  is  not  necessary.  It  is  well  known  (see 
Krein,  page  4 3)  that  one  can  always  introduce  an  equivalent 
norm  so  that  the  factor  M in  (7.4c)  can  be  replaced  by  1, 
however  we  will  eventually  approximate  B by  a family  of 
bounded  operators  B.  where 

h 


(7.5) 


h B,  t 
Zt  E e h 11  i M 


and  we  may  not  be  able  to  introduce  a norm  so  that  M can 
be  taken  as  1 for  all  h,  and  thus  we  will  leave  M unspecified. 

Now  to  make  the  homogeneous  equation  exponentially 
stable  we  introduce  a damping  term  -6(t)  and  consider 
the  equation 


(7.6) 


y = [B  - <5  ( t)  ]y 


For  conditions  on  6 we  first  let  it  be  a scalar  AP  function 


such  that 


(7.7) 


R.P.  (m( 6 ) ) > 0 


This  implies  that  if  g(t,s)  = expj-  J dr  6(r)J  then 


(7.  8) 


g(t,s) I < K e 


-a (t-s ) 


t > s , 
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as  can  easily  be  seen. 


It  is  also  possible  to  permit  6 to  be  an  AP  operator 
function  satisfying  certain  commutivity  properties.  Speci' 
fically  if  Y(t)  is  the  fundamental  solution  to 


(7.9) 

Y = - 6Y 

then 

(7.10) 

(a) 

HY(t)Y_1(s)ll  £ K e"a(t_s)  , t >_  s 

(b) 

6(t),  Y ( t ) , Y ^(t)  commute  with  Zr 

(c) 

for  all  t and  all  r > 0 
6(t),  Y(t),  Y "^(t)  map  Dg  into  itself 

These 

and  commute  with  B. 

hypotheses  permit  us  to  essentially  treat  6 as 

scalar 

in  the 

following  analysis.  If  we  define 

(7.11) 

g( t,s)  = Y ( t) Y~ 1 (s ) 

then  the  "fundamental  solution"  to  (7.6)  is 


(7.12) 


g(t,s)Zt_s  , 


t > s, 


and 


(7.13)  «g(t,s)Zt_s«  <KMe"a(t"S)  EKe“a(t_S) 


after  a redefinition  of  K. 

Now  consider  the  inhomogeneous  equation 


(7.14) 


y = [3-6  ] y + f(t) 


where  f e A(E)  or  C (E) . We  want  to  study  bounded  (for  all  t) 
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1^-1. - 


solutions  to  (7.14).  First  we  see  that  any  such  solution 
is  unique.  In  fact  if  w{t)  were  a bounded  solution  to 
(7.6)  we  would  have  (for  t > s) 


(7.15)  w(t)  = Z g(t;s)  w(s) 

L.“  S 


and  w = 0 follows  immediately  on  letting  s -*■ 

It  is  equally  simple  to  see  that  any  solution  must 
be  given  by 


(7.16)  y(t) 


g(t,s)  z 


t-s 


f(s)  ds 


C(f)  . 


This  is  in  fact  an  immediate  consequence  of  letting  tj. 

in  the  variation  of  constants  formula 

t 

(7.17)  y ( t)  = g(t,tQ) Zt_t  y (tQ)+  j ds  g ( t , s ) Zfc_s f (s ) 

° trt 


which  is  proved  in  Krein  (Theorem  6.1,  page  129)  for  the 
autonomous  case,  and  it  is  a trivial  matter  to  see  that  the 
same  proof  will  work  if  5 is  time  dependent.  Note  that  (7.17) 
need  not  be  a solution  to  (7.14),  but  any  such  solution 
must  be  given  by  (7.17). 

We  now  consider  the  operator  C defined  by  (7.16). 

C is  a bounded  operator  mapping  -*■  . Boundedness  is 

an  immediate  consequence  of  (7.13)  and  in  fact 

(7.  18)  llyll  < K (=  -)  II  fll 

-*  oo  — Cfc  00 


The  integrand  is  obviously  continuous  if  f e and  the 
only  point  in  question  is  its  measurability  for  a general 


f £ L”  (£°r  measurabllity  °f  functions  in  a Banach  space 
see  Hille  and  Phillips  [15]).  However  it  is  shown  by  Krein. 
and  will  be  used  crucially  very  shortly,  that  Zj.  is  the 
strong  limit,  pointwise,  of  where  Bn  is  I bounded 

operator  (Theorem  2 . 9 , page  4 8)  . since  eV  win  be 
norm  continuous  this  establishes  the  integrand  as  the 
pointwise  limit  of  measurable  functions,  hence  measurable. 

We  will  now  show  that  y as  defined  in  (7.16)  is  a 
solution  to  (7.14).  We  will  mahe  use  of  the  above  mentioned 
approximation  theorem  of  Krein  (Theorem  2.9).  It  is  shown 

there  exists  a squence  of  bounded  operators  B^  such  that 

(7.19)  (a)  I (B-Bn)el  — ^ > 0 (»e  s db) 

(b)  I (Zr-2J)e«  0 ve  , r > 0,  where  = SV 

(c)  B Z^ll  <.  m , r > 0 

(d)  Bn  i Zr  (r  0)  satisfy  the  same  commutivity 

relations  with  5(t)  (7.10b,c)  as  does  B and  Z 

r 

We  only  have  to  point  out  that  (7.19d)  follows  from  the 
explicit  definition  of  Bn  as  -AnI  - A^R^  (B  ) where  An  - « 

(see  Krein,  page  49).  We  point  out  that  the  operators  Bfi 
cannot  in  general  be  considered  as  spatial  discretizations. 

A family  of  such  approximations  will  be  assumed  in  the 
following  section. 

Now  consider  the  operators  Cn(f)  defined  by 
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(7.20) 


g(t,s)  Z^_s 


y 


n 


! 


f (s)  ds  EC  (f)  . 

n 


It  follows  immediately  that  the  operators  Cn  are  uniformly 
bounded  operators  mapping  into  and  in  fact 


(7.21)  II  y II  < k II  fll 

n 00  — oo 

where  we  take  the  same  bound  as  in  (7.18).  It  is  also 
immediate  that  yn  is  the  unique  solution  to 

(7.22)  y = [b  — 6 ] y + f 

n n Jn 

and  yn  is  AP  if  f is. 

We  are  going  to  show  that  C -*•  C strongly  in  a certain 
subspace  of  L^.  Specifically  define  S to  be  the  subspace 
spanned  by  the  functions  whose  range  has  compact  closure.  It 
is  well  known  (see  Amerio  [16])  that  ACS.  For  f e s we 
will  have 


(7.23)  Cn(f)  + C(f) 

To  show  (7.2  3)  note  that  by  the  uniform  boundedness  of  C 

n 

it  is  sufficient  to  show  it  for  a dense  subset  of  S and  by 
the  definition  of  S it  is  sufficient  to  take  f as 

(7.24)  f = l x (t)e  . 

j=l  3 3 

where  e . e e and  Xj  is  the  characteristic  function  of  a 
measurable  set.  By  linearity  we  only  have  to  work  with 
one  term  in  the  sum  in  (7.24),  Thus  if  I is  some  measurable 
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set  let  f - Xle.  Then  wn(t)  = y(t)  - yn(t)  is  explicitly 


(7.25)  wn(t)  = J x(t-r)  g(t,t-r)  [Zr-Z^] 


e dr  . 


Estimating  (7.25)  we  obtain 


00 

(7,26)  11  wn  11  oo  - K / e~ar|1  (Zr-Z^)ell  dr 

0 

and  (7.23)  is  an  immediate  consequence  of  the  Lebesgue 
dominated  convergence  theorem. 

Observe  that  (7.23)  implies  that  C(f)  is  continuous 
if  f s s and  it  is  AP  if  f is. 

Now  consider  the  almost  periodic  case.  First  we  point 
out  that  if  6 and  f are  C*  (AP  derivatives)  then  y is 
also  and  in  fact 


(7.27) 


y = C (-6y  + f) 


This  is  trivial  if  we  observe  that 


(7.28) 


yn  " Cn'-6yn  + 


and  if  y is  defined  to  be  the  right-hand  side  of  (7.27)  we 


have 


(7.29)  y-Yn-  IG-Cn]  <-Sy)  + Cn(-6  (y-yn)  ) + [C-Cn]  (f)  , 


and  clearly 


(7.30) 


11  (y-yn)  II  > 0 

2 n oo  n-K>o 


Equation  (7.27)  is  now  an  immediate  consequence  of  lettii 
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n + “ in  the  equation 


(7. 31) 


Yn  (tl)  ~ Y-  (t'>)  = Y-  (t)  dt 


n 2 


n 


We  next  show  that  y satisfies  an  integrated  version 
of  (7.14).  Specifically 


(7.32)  y(t1)  - y(t2 


r r f2 

= B j dhy(t)  — J dt  6 ( t)  y ( t ) + I dt  f ( t) 


Since  is  dense  and  B is  closed  it  is  sufficient  to  show 

(7.32)  when  f is  of  the  form 


(7.33)  f = b , be  dd. 

B 

Now  observe  that  Bf  is  AP.  It  certainly  follows  that 
(7.  34)  By  = C (Bf ) . 


In  fact  this  holds  whenever  Bf  e CM  and  is  a simple  conse- 
quence of  the  closure  of  B,  approximating  the  integral  in 
(7.16)  by  Riemann  sums  over  finite  intervals  and  then 
noting  that  B applied  to  each  sum  is  the  Riemann  sum 
approximating  C(Bf)  . 

Now  the  same  property  certainly  holds  for  B and  C i.e. 

1 n n 


(7.35) 


B y 
nJ  n 


= C (B  f ) . 
n n 


Finally  since  by  its  construction  Bn 


commutes  with  Z 

r 
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(Krein,  Theorem  2.9)  we  can  write 


(7.  36)  B y = C (B  f)  . 

nJ  n 

Since  yn  satisfies  (7.22)  we  have 


(7.37)  yn(tl)  - yn(t2) 


j 2 ^2  ^2 

= Bn  J dt  yn(t>  - f dt  6(t)yn(t)+  J dt  f(t) 


= B 


[2 

t0  t 

f2 

dt  y ( t)  + [B  -B] 

J n 

dt  y(t)  + B 
J n J 

dt  (yn-y) 


! )2 

dt  6 (t)yn (t)  +|  dt  f (t)  . 


where  the  fact  that  j dt  y(t)  £ D0  is  a consequence  of 
(7.34)  and  the  same  type  of  closure  argument  as  followed 
that  equation. 

If  we  now  let  n ->■  °°  in  the  last  equation  in  (7.37)  we 

[*~2 

see  that  the  term  Bn  I dt  (yn~y)  — — > 0 because 


(7.38)  Bn(yn-y)  = [Cn“C](Bnf)  = [Cn~C](Bf)  + (C^C ] ( (Bn~B ) f ) 


and  the  right-hand  side  of  (7.38)  clearly  0 in  the  sup  norm 
as  n -*  °°.  It  therefore  follows,  on  letting  n “ in  (7.37), 
that  (7.32)  does  in  fact  hold  for  functions  of  the  form  (7.33) 
and  thus  for  all  f e A. 

Now  if  6 and  f are  C (hence  y is)  we  can  show  very  easily 


that  (7.32)  implies  (7.14).  In  fact 
let  h 0 in  the  equation 


y ( t+h)  -y  ( t) 
(7.39)  r 


t+h 

B | ds  y(s) 
t 

h 


since  y is  C we  can 


t+h  t+h 

ds  6 ( s ) y (s ) ( f ( s ) 


where  h can  be  positive  or  negative,  using  the  closure  of  B. 
Similarly  (7.14)  will  hold  if  we  only  require  Bf  is  AP. 

In  fact  since  in  that  case  By  is  AP,  B can  be  brought  inside 
the  integral  in  (7.32)  , which  can  then  be  differentiated 
directly. 

We  point  out  that  in  both  these  cases  both  y and  By  are 

also  AP.  For  future  use  we  point  out,  using  the  integral 

2 

equation  (7.32)  applied  to  (7.27),  that  if  6 and  f are  C 
2 

then  y is  C and  y satisfies 

• • 

(7.40)  y = [ B— 6 ] y - 6y  + f 

i.e.  (7.14)  can  be  differentiated  formally.  (Note  in  parti- 
. 2 

cular  that  if  Bf  is  AP , B y is  AP.)  This  argument  can 
obviously  be  extended;  namely  if  f and  6 are  then  y is 
and  we  have 

(7.41)  y r = [B-  <5  ] y r 1 + fr  1 - (6y)r  1 + 6yr  1 , r=l,...,p. 

We  will  now  examine  briefly  the  case  that  f e S. 

• • 

We  still  have  y -*•  y but  (7.30)  need  not  hold  since  6y  and  f 
need  not  belong  to  S.  However  we  see  very  easily  that  if 

f e L and 

00 
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(7. 42) 


x = C(g) 


Xn  = Cn(g} 

then  we  will  have  xr  ( t)  - x(t)  pointwise  and  boundedly.  In 
fact  if  wn(t)  = x(t)  - x (t)  , then 

oo  * 1 

(7.43)  lwn(t)l  < K I e'“rl(Zr-ZJ)  g (t-r)  II  dr 

0 

and  the  Lebesgue  dominated  convergence  theorem  is  still 
applicable  and  hence  (7.31)  will  still  yield  ( 7.27 ),  if  we 
use  the  bounded  convergence  theorem  and  the  fact  that 

(7.44)  y(t)-yn(t)  = fOC  ] (-6y+f)  (t)+  C (-6(y-y  ))(t)  > 0 

n n n-*00 

Thus  y has  a bounded  derivative.  Now  if  we  assume  that  f is 
of  the  form 

(7.45)  f = y b > 

r XID  ' b e Dfi  , 

then  the  derivation  of  (7.32)  from  (7.37)  and  (7.38)  still 
follows,  if  we  use  the  bounded  convergence  theorem  for  inte- 
grals over  finite  intervals.  Since  finite  sums  of  such  f's 
are  dense  m S,  (7.32)  holds  for  all  f e s and  if  Bf  e S or 
A and  f are  C1,  y will  still  be  the  unique  Cm  solution  to 
(7.14).  Thus  the  theory  will  hold  for  a certain  class  of 
bounded,  non  AP , inhomogeneous  terms,  but  is  not  as  simple 

as  the  AP  case  and  in  the  future  we  will  restrict  ourselves 
to  this  case. 

Finally  we  point  out  that  (7.32)  can  be  derived  by 
brute  force  dif ferentation  of  the  formula  (7.16)  but  use  of 
the  approximating  operators  Bn  permits  a cleaner  and  more 
straightforward  development. 
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8*  Discretization  in  Space 


The  operators  which  proved  so  useful  in  the 

preceding  section  need  not  correspond  to  a discretization 
in  space.  We  therefore  postulate  the  existence  of  a 
family  of  bounded  operators  Bh  , where  h will  generally 
play  the  role  of  a spatial  grid  size,  such  that: 


(8.1)  (a) 

(b) 

(c) 


Bh  are  ^efi-ned  and  bounded  for  h e (O/hg] 

"V  = ' * > ° 

h 

Bh  generate  bounded  semi-groups  (i.e. 
Bht  . 

He  = Zall  < M for  t > 0) 


(d)  e G D implies  II  (B-B.  ) ell  — — — ->  o 

B h h+0 

(e)  If  6 is  not  a scalar  function,  then  the 

commutivity  relations  (7.10)b,c  hold  with 

B and  z replaced  by  B,  and  Zh 
r hr 

(f)  If  e e (or  any  dense  manifold  in  Dfi 
invariant  under  Z^  for  r >_  0)  we  have 

llBhZrel1  - K(I'e^  for  r 6 I,  where  I is  any 
compact  interval  on  the  nonnegative  real  axis. 


We  point  out  that  (8. If)  is  designed  to  insure  that 
Zr  Zr  stron9ly#  for  all  r > 0.  This  is  certainly  reasonable 
since  B^Z^e  converges  to  BZ^e.  it  is  of  course  trivial  if 
Bj-j  commutes  with  Also  note  that  the  convergence  in  (8.  Id) 

does  not  deal  with  the  order  of  the  approximation  of  B,  to  B. 


■e—.Ui-.jIi.i.  1,.,,,..  ; ..  


This  will  be  dealt  with  by  a subsequent  assumption. 

Before  proceeding  further,  we  indicate  a simple  example 
of  the  case  we  are  considering. 

(8.2)  (a)  E = L2  ( 0 , 2tt  ) 

(b)  b = — (with  periodic  boundary  conditions) 

(c)  B = u ( x+h) -u(x-h) 

h 2h 


In  this  example  (8. If)  is  unnecessary  as  Bh  in  fact  commutes 

wi th  Z . 
r 

Now  associated  with  B^  we  have  the  equations 


(8-3)  y = [Bh-6]y 

(8-4)  y = [Bh-fi]y  + f 

The  homogeneous  equation  (8.4)  has  as  its  fundamental  solution 

Bhfc 

'* h ( = g(t's)  e , and  is  exponentially  stable,  i.e.  there 
exist  constants  K and  a independent  of  h such  that 


(8.5) 


ilYh(t)Y^1(s)ll  < K e"a(t_s)  , t > s. 


Associated  with  (8.4)  we  form  the  operator  where 

t 

(8.6) 


yh  = ch(f)  = J g(fc's)  zt-s  f(s)  ds 


for  any  AP  f,  yh  is  the  unique  AP  solution  to  (8.4). 
Equation  (8.5)  shows  that  the  operators  are  bounded 

uniformly  in  h , i.e. 


(8.  7) 

where  K is  independent  of  h 


yjl  < K II  fll 

n o°  — oo 
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Now  our  first  task  is  to  show  that  Ch  ■+•  C strongly. 

Using  the  results  of  Section  7 (see  (7.23)  ff)  it  is  suffi- 

h 

cient  to  show  that  Z^_  -*•  Z^_  strongly,  for  r > 0.  This  will 

be  an  easy  consequence  of  (8. If).  In  fact  if  e £ d (or 

B 

the  manifold  described  in  (8. If))  and  if  y = Z e , 
h h 

Yh  = zte  / then  if  w = y-y^  we  will  have 


(8.8) 

h „ 

w = B,  w, 
h h 

h 

+ [B-Bh]y 

wn(0)  = 

0 . 

From  (8.8) 

we  can  write 

t 

, f 

(8.9) 

w (t)  = J ds 

Zts[B-V2se 

i h 

and  the  fact  that  II  w (t)U  ->0,  in  fact  uniformly  on  compact 

intervals,  is  an  immediate  consequence  of  (8. If)  and  the 

Lebesgue  bounded  convergence  theorem.  The  result  for  any 

e € e follows  from  the  uniform  boundedness  of  the 
h - 

operators  Zr  and  the  denseness  of  the  manifold  of  (8. If). 
Wow  the  fact  that 


<e.).(»  "<y-yh>"»  -h^->  0 

for  all  f £ A (in  fact  for  all  f £ S)  is  not  any  improvement 
of  the  corresponding  result  for  the  operators  Bn  of  Section  7. 
We  want  to  show  that  the  difference  in  (8.10)  approaches  zero 
as  fast  as  some  power  of  h,  and  to  do  this  we  must  require 
that  B^  approximate  B up  to  a certain  degree  of  accuracy. 
Specifically  we  suppose  that  there  exists  an  operator  L 
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(unbounded)  with  D^  c db  and  such  that  if  e t DL  we  will 
have 


(8.11)  H(B-B^)eU  c h^llLell  where  j > 0 

and  C and  j are  independent  of  h and  e. 

If  we  have  Lf(t)  AP  we  would  like  to  show  that  (8.11) 
implies 

(8-12)  II  (y— yh)  II  ^ = 0(h:i)  . 


Unfortunately  (8.12)  requires  a further  assumption;  namely 


(8.13) 


(a) 

(b) 


If  Lf  is  AP  then  IlLZ  f(t)l  will  pe 
uniformly  bounded  for  r > 0 and  all  t. 


Assumption (8.13)  is  of  course  trivial  if  L commutes  with  Z . 

r 

It  is  now  a simple  matter  to  show  (8.12).  in  fact 
wh  = y~Yh  is  the  unique  AP  solution  to 

(8-14)  wh  = [Bh-6]wh  + (B-Bh)y 

whence  by  (8.7)  we  have 


<8-i5)  '^h11®  = 0(11  [B-Bh]yll  J 

and  so  it  is  only  necessary  to  study  (B-Bh)y.  But  now  if 
Bf  e a,  which  will  generally  follow  if  Lf  e A since  L will 
be  a higher  order  operator  than  B,  we  can  write 
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(8.16)  [B-BhJy  = j dr  g ( t , t-r ) [B-Bh 

0 


and  clearly 


(8.17) 


11  [B-Bh]y“co  = °(h3) 


]Zrf (t-r) 


is  an  immediate  consequence  of  (8.13). 


9 ' Discretization  in  Time  by  the  Multistep  Linear  Method 


We  will  now  apply  the  theory  developed  in  Section  3 
to  the  spatially  discretized  equation 

(9-1>  yh  = fBh"5]  + f • 

Since  Bh  is  bounded  the  theory  is  immediately  applicable  and 
we  can  assert  that  for  k sufficiently  small  (but  depending 
on  h)  there  will  be  an  AP  sequence  wR(k,h)  which  will  converge 
to  wn(k,h)  defined  by 


(9.2) 


wn  (k  ,h) 


f Yh  ( (n+£-l)  k)  • 
' yh(nk) 


and  if  everything  is  smooth  in  t we  will  have 


(9.3) 


lwn-wn“~  = °h(k  > 


where  the  bound  in  (9.3)  may  depend  on  h.  It  is  only  left 
to  show  that  the  bound  in  (9.3)  is  independent  of  h and  to 
discuss  the  relationship  of  k to  h.  We  will  show  that  what 
is  required  is 


(9- 4)  k + "4r  = 0(1) 

li 

i.e.  there  exists  a constant  c independent  of  k and  h, 
such  that  if  k + k/h  m < c the  difference  equation  corres- 
ponding to  (9.1)  will  be  exponentially  stable  with  constants 
independent  of  k and  h. 
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. • t:  ‘ 


The  quadratic  dependence  on  h in  (9.4)  is  restrictive 
and  will  be  removed  in  the  following  section. 

Throughout  this  section  we  will  use  the  notation  of 
Sections  2 and  3.  We  will  also  take  6 to  be  C^. 

We  now  fix  h and  consider  the  homogeneous  difference 
equation  on  (see  (3.7)), 

(9.5)  w = U(kB  - k6  ) + k20(l)Wri  . 

n+1  h n n 

It  is  clear  that  this  will  hold  uniformly  in  k and  h if 
k + k/hm  is  sufficiently  small,  i.e.  if  we  are  in  a region 

(9.6)  k + = 0(1)  . 

hm 

Of  course  restrictions  of  the  form  (9.6)  are  included  in  (9.4). 
Also  note  from  the  explicit  definition  of  the  full  linear 
term  in  (9.5)  (see  (2.11))  the  last  term  in  (9.5)  will  be 
uniform  in  k and  h.  For  the  rest  of  this  section  this  will 
be  understood  for  all  "0"  signs  unless  stated  otherwise. 

If  Wn(k,h)  is  the  fundamental  solution  to  (9.5)  it  is 
clear  that  all  we  must  show  is  that  in  a region  of  the  form 
(9.4)  we  will  have 

, -a  k(n-j) 

(9.7)  II W W . It  < K.  e 1 , n > j , 

where  the  constants  and  are  independent  of  k and  h. 
Following  the  procedure  of  Section  3 we  let 

(9.8)  v = T(kB,  - k6,  )w 

n h li  n 
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and  we  then  have  in  a 


region  of  the  form  (9.6) 


(9.9) 


Vn+1  = L(kBn  " k6n)Vn  + k 0(1)vn 


Using  the  structure  of  the  block  diagonal  operator  L(z) 
derived  in  Theorem  1 we  can  write  this  as 


(9.10)  v 


n+1 


[I+k  (B  -6  ) 0 

n n 

0 B 


n 


v + kO(k  + -9-)v 
n ^m  n 


Here  we  have  used  the  notation  7 for  the  lower  block  (see 

(2.26))  so  as  not  to  confuse  it  with  the  operator  Bh . 

Observe  that  in  (9.10)  we  have  the  first  appearance  of 
2 rn 

the  term  k/h  . This  is  because  in  order  to  apply  Theorem  2 
we  must  separate  out  a factor  of  k in  the  perturbing  term. 

How  it  certainly  follows  that  in  a region  of  the  type 
(9.6)  we  have 

(9.H)  II  77  B II  < K en'j  , n > j, 

i = j 

where  K1  and  0 are  independent  of  h and  k.  We  therefore 
only  have  to  consider  the  upper  block  in  the  leading  term 
of  (9.10)  and  in  order  to  apply  Theorem  2 we  are  going  to 

compare  this  with  Yn;h  = Yh(nk)  where  Yh  is  the  fundamental 
solution  to 


(9-12)  *h  - iVslYh 

Observe  that  the  estimate 

(9.13)  II Y ,Y-1  , H < K e-ak(n-j) 

n ; h 3 ; h — • 
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holds,  where  K and  a are  independent  of  h. 


We  will  show  that  Y , satisfies  the  difference 

n;h 


eq uation 


(9.14) 


Y . . , = [I+kD  (k,h) ] Y , 

n+1 ;h  n n ; h 


= [ I+k  (B.  -6  ) J Y , + kO  (k  + -~)Y  , . 

h n n;h  h2m  n;h 


Observe  that  the  last  term  in  (9.14)  is  of  the  same  order 
as  the  perturbation  in  (9.10). 

Equation  (9.14)  follows  easily  from  the  same  argument 
as  in  Section  3 (see  (3.13)  and  (3.14))  if  we  simply  observe 
that 


(9.15) 


(Bh  -6)  = - 6 


and  that  (3.15)  holds  in  this  case  for  t > s by  (8.5). 
We  can  now  apply  Theorem  2 comparing  the  equation 


(9.16) 


V*  , = [I+k(B  -5  )]V* 

n+ 1 ; h h n n ; h 


with  (9.14).  We  simply  note  that  the  proof  of  Theorem  2 


indicates  that  the  terms  and  (see  the  statement 


of  Theorem  2)  depend  only  on  the  constants  K and  u in  (9.13) 
and  are  thus  independent  of  both  k and  h.  Another  applica- 
tion of  Theorem  2 comparing  the  equation 


(9.17) 


V 


r I+k<VV  0 


n+1 


B 


V 


n 


n 


with  (9.10),  together  with  the  boundedness  of  the  operators 
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T(k(Bh~6n))f  T ^(k(B^-6  ))  in  the  region  described  by  (9.6) 
will  yield  (9.7)  with  the  restriction  k + k/h2m  sufficiently 
small  (i.e.  (9.4)). 

If  we  now  refer  to  the  proof  of  admissibility  in  the 
stable  case  given  by  (3.31)  ff  we  see  immediately  that  (9.7) 
implies  (L^,!^)  admissibility  for  the  inhomogeneous  version 
(9.5),  with  an  admissibility  bound  independent  of  k and  h. 

It  is  only  necessary  to  show  that  the  convergence  of  w^  to  w^ 
is  independent  of  h (see  (9.3)).  This  will  follow  from  (3.40) 
if  we  can  show  By^+^!loo  will  be  bounded  uniformly  in  h, 
assuming  that  6 and  f are  smooth. 

In  fact  if  6 and  f have  p+1  derivatives  the  differenti- 
ability properties  of  y^  are  the  same  as  those  of  yn  introduced 
in  Section  7.  In  addition,  the  formula  given  in  (7.28)  can 
obviously  be  extended  (for  the  operators  C^)  to  yield 

(9.18)  yP+1  = Ch(fP+1  + (-6yh)P+1+  6yP+1) 


(note  that  the  argument  of  in  (9.18)  does  not  involve  yp+^) 

and  the  uniform  boundedness  of  the  operators  immediately 

yields  the  uniform  boundedness  of  the  derivatives  of  y . 

r y ( (n+&-l)k  - h 


If  we  now  define  w = 

n 


l 


, then  (9.3) 


y (nk) 

(uniformly  in  h)  together  with  (8.12)  implies 


(9.19)  **wn-wn**°°  = °(kP+h^)  = 0(h2mp+h^) 

where  we  have  used  (9.4)  for  small  k. 
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10 . The  Lax-Wendroff  Scheme 
2m 

The  term  k/h  in  (9.4)  can  be  improved  for  operators 
such  that  the  Lax-Wendroff  scheme  is  stable  for  the  homogeneous 
equation 

(10.1)  w = Bw 

We  will  then  be  able  to  obtain  AP  sequence  solutions  under 
the  more  favorable  restriction 

(10.2)  k + ^ = 0(1)  . 

First  of  all  we  require  that  6 and  f are  C3.  As  shown 
in  Section  7 this  implies  that  y is  C . It  implies  further 
that  y satisfies 

(10.3)  y = [B-6]y  - 5y  + f 
i.e.  that  the  basic  equation 

(10. 4)  y = [b-6  ] y + f 

can  be  differentiated  formally. 

We  can  now  apply  the  Lax-Wendroff  scheme  to  (10.4). 
Proceeding  formally  we  write,  assuming  Bf  is  AP, 

k2 

(10.5)  y ( t+k)  = y { t ) + ky  ( t)  + — y(t) 

k2  2 

= y (t)  + k [B-6  ) y ( t)  + ~ [B-6]  y(t) 
k2  ' 

+ — [-6y  + (B-6 ) f +f]  + kf  ( t)  = 
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- [I  + kB  + y-  B‘']y(t)  - k6y 
+ k[-y-  62-  | 6 - k6B]y+  k[f+  ~(B-6)f+  | f] 


We  have  written  this  as  an  inhomogeneous  difference  equation 
in  order  to  apply  techniques  already  developed. 

Now  in  the  expression  for  y we  approximate  B by  B 

h 

and  in  the  expression  for  y approximate  B2  by  B2  , for  a 

ch 

certain  constant  c > 0 to  be  determined  (c  = 1/2  for  the 

example  described  in  (8.2)).  Of  course  we  will  have  to 
2 o 

assume  that  B.  approximates  B in  the  same  sense  that  B 

h 

approximates  B,  but  we  will  leave  the  precise  hypothesis 
for  later. 

Applying  these  approximations  to  (10.5)  and  letting  w(k,h) 
stand  for  the  dependent  variable  we  obtain  the  difference 
equation 


(10.6)  wn+1(k,h)  = u(k,h)wn-  k6nwn+  kR(n,k,h)wn+  kf 
where 

.2  - 

(10.7)  (a)  U(k,h)  = I + kB,  + B^u 

h 2 ch 

(b>  'fn=  f„  + ! lBch  - + I K 

(O  R(n.k.h)  -|^-|«n-k6nB ch  = 0(k+^). 

h 

Now  U(k,h)  represents  an  approximation  to  (10.1).  We 
choose  c so  that  for 

k 


(10.8) 


m 


= 0(1) 
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this  discretization  is  stable,  i.e. 


(10.9)  I U(k,h) rl  < H , r > 0. 

We  are  going  to  use  Theorem  2 to  show  that  if  HR(n,k,h)l 
(i.e.  k + k/hm)  is  sufficiently  small,  then  the  homogeneous 
version  of  (10.6)  is  exponentially  stable;  that  is  if 
Wn(k,h)  is  the  fundamental  solution  to  the  homogeneous 
equation,  then  there  exists  constants  and  independent 

of  k and  h such  that 

. -ot.k  (n- j ) 

(10.10)  «wnWj  I £ K1  e , n > j, 

where  for  simplicity  we  have  suppressed  the  dependence  of 

on  k and  h.  Once  we  have  obtained  (10.10),  (L^,!^)  (and 

(A  ,A  ) ) admissibility  will  follow  in  the  usual  manner, 
n n 

To  obtain  (10.10)  we  observe  that  since  we  can  allow 
k + k/hm  to  be  as  small  as  required,  Theorem  2 permits  us 
to  neglect  the  perturbing  term  R.  Also  by  neglecting  a 
perturbation  of  the  same  order  as  kR  we  can  work  with  the 
simpler  equation 

(10.11)  w^+1  = [I-k6n)  U(kh)w*  . 

Now  consider  the  function  g(t,s)  defined  in  (7.11). 
Using  the  same  argument  given  in  Section  3 (see  (3.16)),  we 
see  that  gn  = g(nk,0)  will  satisfy  the  equation 

(10.12)  gn+1  = [I+kDn]gn  = [I-k6n+k20 (1) ]gn 

where  the  "0"  sign  in  (10.12)  is  independent  of  n.  Thus  by 
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aqain  neglecting  a perturbation  of  the  form  kO(k+k/hm) 
we  obtain  finally  the  equation 


(10.13) 


wn+i  = U+kDn]  U(k,h)wn  . 


But  the  fundamental  solution  to  (10.13)  is  simply 


(10.14) 


-1 


wnwj  = g(nk,  jk)u(k,h)n_:) 


as  one  can  verify  immediately,  using  the  commutivity 
properties  of  gn  with  Bh  (see  (8.1e))  in  the  case  that 
6 is  an  operator.  Now  the  estimate 


(10.15) 


_ _ -ak(n-j) 

I W W . I < i<  e 
n j — 


n 1 3 


follows  from  (7.10a)  and  (10.9)  and  we  can  now  apply 
Theorem  2,  noting,  again  that  the  terms  K^,  a^,  depend 

only  on  the  terms  K and  a (see  the  hypothesis  of  Theorem  2) , 
to  obtain  (10.10) . 

We  now  have  (1^,1^)  (or  (A^A^))  admissibility  for 
(10.6)  and  it  is  only  necessary  to  study  the  convergence  in 
the  case  that  the  inhomogeneous  term  is  given  by  (10.7b). 

To  do  this  we  must  now  assume  that  B2  approximates  B2 ; namely 
that  there  exists  an  operator  L such  that  if  e S D n d 


then 


B 


(10. 16) 


UB2-B2)ell  < C2hj#£elf 


We  also  assume  that  L satisfies  the  same  assumption  as  L 
(see  (8.13))  so  that  if  Lf  and  Lf  are  AP  we  can  conclude 
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(using  the  same  argument  as  in  (8.16)  and  (8.17)) 


(10.17)  (a)  I (B-B,  )yl  = O(h^) 

h J *> 

(b)  II  (B2-B2)y#  x = 0(hD)  . 

If  wn(k,h)  is  the  unique  AP  solution  to  (10.6),  we  can  now 
show  that  w^  converges  to  yn  uniformly  in  n. 

3 

In  fact  since  y is  C , yn  satisfies  (10.5)  up  to  an  error 
3 

0(k  ).  If  we  then  make  use  of  (10.17)  we  see  that  y will 

Jn 

satisfy  (10.6)  with  an  error  kO(h^+k2),  where  this  is  uniform 
in  n.  Then  as  the  admissibility  bound  of  (10.6)  is  independent 
of  h (by  (10.10)),  we  can  conclude 

(10.18)  “wn-ynloo  = 0(hj+k2)  = 0(hj+h2m) 

where  we  have  used  (10.2). 

For  the  example  described  in  (8.2)  we  will  have 

33  - a4  . 

(10.19)  L = — ^ and  L = — j , c = ^ , m = 1,  j = 2 , 

3x  9x 

(periodic  boundary  conditions) 
Finally  we  would  like  to  point  out  that  in  the  term  f 

n 

in  (10.7b)  one  can  replace  Bchfn  without  changing 

the  order  of  the  error.  This  is  simply  because  Lf  is  AP, 
hence 

(10.20)  II  (B-B.  ) f H = 0(hj)  . 

h °° 
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11.  Linear  Perturbations 


We  now  consider  the  perturbed  system 


(11.1) 


y = [B-<5  ] y + D ( t)  y + f(t) 


Here  D(t)  is  a bounded  operator  for  each  t and  is  almost 
periodic  as  an  operator  function.  We  also  require 


(11.2) 


Di  < 

OT  — 0 


where  Cq  is  to  be  specified. 


The  first  step  in  finding  AP  solutions  to  (11.1)  is 
to  solve  the  equation 


(11.3) 


y = C (Dy ) + C(f) 


This  equation  can  be  solved  uniquely  by  the  contracting  map- 
ping principle  provided  u = GqK  < 1,  where  K is  the  norm 
of  C (see  (7.18)).  The  contracting  mapping  principle  also 
shows  that  the  assignment  of  y from  f is  a bounded  operator 
and  in  fact 


(11.4) 


y II  m II  fl  = k,I  fll 

■<  00  _ 1-  y 00  1 


where  depends  only  on  and  is  independent  of  D. 

We  want  to  show  that  y satisfies  (11.1),  and  if  we 
assume  that  D,  6 and  f are  C1,  the  results  of  Section  7 
show  that  it  is  sufficient  to  show  that  y is  C^.  This  will 
follow  if  we  observe  that  y will  be  the  limit  of  iterates 
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1 


0 


defined  as 

(a) 

(11.5) 

(b)  yi+1  = C(Dyi)  + C(f)  . 

We  will  then  have  for  each  i 

(11.6)  y.+1  = C(Dy.)  + C(-6Yi  + Dy.  + f) 

If  we  then  define  y as  the  unique  solution  to 

(11.7)  y = C (Dy ) + C(-6y  + Dy  + f) 
we  will  have 

(11.8)  yL  - y. 

To  see  (11.8)  note  that  H y± H ^ is  bounded  uniformly 
in  i (from  (11.6)  and  u < 1)  and  if  g^^  is  defined  as 

(11.9)  = sup  H y-yr» 

r>i 

Then  {g^}  form  a nonincreasing  sequence  satisfying 

(11.10)  g.+1  < u + o(l) 

^403 

and  (11.10)  implies  g^^  -*•  0 which,  in  turn,  implies  (11.8). 
The  fact  that  y = y follows  exactly  as  in  Section  7 (see 
(7.27)  ff)  . Observe  that  since  y is  AP,  By  is  also  AP 
( from  (11.1)). 

Wow  consider  the  operators  Ch  defined  by  (8.6).  The 
equation 

(11.11)  yh  = ch(Dvh>  + ch(f> 


-97- 


can  be  solved  exactly  as  (11  3)  can.  (For  simplicity  we 

suppose  that  Icy  = ic«  sc  that  we  have  the  same  restrio- 
tion  on  eQ.) 

Our  first  task  is  to  show  that  y^  -*  v.  This  will  be  a 

simple  consequence  of  the  strong  convergence  of  C to  C 

h 

In  fact  if  w^  = y - we  have 

(11.12)  wh  = (C-Ch](Dy)  + ch(°wh)  + [C-Ch ] ( f ) 
whence 


(11.13) 

and  this  certainly  implies  II  w 


1 Wh“  00  1 + °(1) 

h-*-0 


h-*0 


->  0, 


Wow  of  course  this  would  not  be  useful  unless  we  could 


prove 


(11. 14) 


lw,  I = 0(hD) 


This  estimate  can  be  shown  provided  we  make  certain 

assumptions  on  the  interaction  of  D and  the  operator  L 

defined  in  (8.11).  We  first  observe  that  w,  satisfies 

h 

the  equation 


(11.15) 


wh  = ch<Dwh>  + ch(IB-Bhly) 


and  using  the  analogue  of  (11.4)  for  the  equation  (11.11), 
together  with  the  uniform  boundedness  of  the  operators  C 


we  see 


(11.16) 


Voo  = 00  [B-Bjylj  , 
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and  it  is  sufficient  to  show 

(11.17)  I [B-B,  ]yl  = O(h^)  . 

n °° 

To  obtain  (11.17)  we  make  an  assumption  on  the  operator 
L.  We  will  call  it  Assumption  I as  we  will  later  have  to  make 
the  same  type  of  assumption  for  other  operators. 

Assumption  I:  Suppose  there  exists  a sequence  of 

operators  LQ  = I,  1^ = L such  that 


(11.18)  (a) 

(b) 

(c) 


(d) 


are  closed 

E = Dt  3D.  3 ...  3 DT 

ii/>  i-i  i — L 

0 1 n 

L.  commutes  with  Z for  t > 0 (i.e.  Z,  : DT  ->DT 
^ t t Li . L . 

1 1 

and  L.Z  e = Z L.e  for  e e DT  ) 

1 t t 1 Li  , 

1 

If  e e D then  D(t)e  e D_  and 
^ 0 


L D(t)e  = l E*(t)L.e 
* j=0  3 J 

o 

where  the  operators  E.(t)  are  bounded  and  AP, 
o 

and  E^  = D. 


Assumption  I is  clearly  motivated  by  the  example  where 
D ( t ) is  multiplication  by  some  function  d(t,x)  (x  is  the 
spatial  variable)  . 

Now  Assumption  I implies  that  if  L^f  is  AP  then  L^y  is  AP, 
for  £ = l,...,n.  First  note  that  if  z = C(f)  then 


(11.19)  L^z  = C(L£f)  . 

This  is  a consequence  of  (11. 18a, c)  using  the  same  argument 

as  was  used  for  B (see  (7.34)  ff).  Now  if  we  define  z„  .=  L„y. 

x.;i  i 

where  the  yi  are  defined  in  (11.5);  then  it  certainly  follows 
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from  (11.18d)  that  is  AP  for  i and  for  i = 0,...,n. 

We  will  now  show  that  the  sequence  {z„.}  converges  as  i->  *>. 

X/  ; 1 

To  see  this  first  set  l = 1 . We  then  have 


(11.20) 


‘l ; i+1  = C(Dz1;.)  + C (E  q ( t ) y i ) + C(f) 


and  zl,i  Z1  ' where  zi  is  the  solution  to 
(11-21)  z1  = C(Dz1)  + C(Ejy)  + C(f)  . 

In  fact  wi  i = zi  “ i satisfies  the  equation 

(11.21)  wi,i+i  = c(Dwi,i)  + C (Eq ( y — y i ) ) 

and 


(11.22) 


1 w . II 

l,i  00 


00 


> 

2_  -kx> 


0 


follows  from  the  same  argument  used  in  proving  y is  C1 
(see  (11.9)  ff ) . Finally  by  the  closure  of  we  must 
have  z ^ = L^y.  It  is  clear  that  a simple  induction,  using 
the  same  proof,  will  establish  that 


(11-23)  Vi  ~T^r>  V ' Z = 1 n- 

In  particular  L (=  Ln)y  is  AP  and  (11.17),  hence  (11.14) 
is  valid. 

Finally  we  point  out  that  i f we  assume  a chain  of 
operators,  as  in  Assumption  I,  for  the  operator  B,  we  will 
then  have  BDy  AP.  We  will  also  have  to  make  such  an 
assumption  for  the  operator  L (see  (10.16)),  when  applying 
the  Lax-Wendroff  scheme  to  (11.1),  as  will  be  discussed 
shortly . 
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We  will  now  apply  to  (11.1)  the  two  methods  of  time 
discretization  that  have  already  been  described.  We  are 
interested  in  the  existence  of  AP  sequence  solutions,  or 
equivalently  exponential  stability  for  the  homogeneous 
difference  equation,  with  k and  h restricted  according  to 
relations  (9.4)  and  (10.2).  As  one  might  expect  this  will 
be  a simple  consequence  of  Theorem  2,  provided  Gq  is 
sufficiently  small. 

We  consider  first  the  linear  multistep  scheme  as 

discussed  in  Section  9.  Since  the  linear  term  C of  the 

n 

difference  scheme  (see  (2.11))  is  a smooth  function  of  its 
£+1  arguments,  it  follows  that  if  k + k/hm  is  sufficiently 
small  (i.e.  IkB^I  is  small)  then  the  difference  equation 
for  (11.1)  can  be  written 

(11.24)  w . = C w + k 0 ( c n ) w + k f 

n+1  n n On  n 

* 

where  is  the  linear  term  corresponding  to  the  unperturbed 

version  of  (11.1),  i.e. 
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Of  k and  h (see  (9.7)  ff)  it  will  follow  from  Theorem  2 

(see  the  remark  following  (9.16))  that  in  the  same  region 

the  homogeneous  version  of  (11.24)  will  be  exponentially 

stable  if  e0  is  small  enough.  It  then  follows  easily 

that  the  inhomogeneous  equation  (11.24)  will  have  (Ln  Ln) 

~n  ~n  00 ' oo ' 

(or  (A  ,A  ))  admissible,  with  an  admissibility  bound 

independent  of  k and  h. 

We  can  thus  obtain  AP  sequence  solutions  to  (11.24) 
and,  as  shown  in  Section  9 (see  (9.18)  ff) , we  will  obtain 
convergence  of  Oik?)  provided  lyP^  is  bounded  in  h. 
However,  as  one  can  trivially  see,  if  D,  6,  and  f are  Cp+1 
then  the  same  argument  used  to  show  y is  (see  (11.6)  ff) 
can  be  repeated  p+1  times  and  applied  to  the  operators  C 

h 

as  well  as  C to  show  that  y and  yh  are  CP+1  and  in  fact  (11.7) 
can  be  generalized  to  (compare  with  (9.18)) 

(11.28)  y1  = C (Dy1)  + ClKD-jJy]1  - (D  - 6)yi  + f1)  , 

1 1 / • • • / P+ 1 f 

with  a similar  formula  for  yh  and  Ch*  This  certainly  shows 

that  Ily^D  will  be  bounded  in  h for  i = l,...,p+l  and 

thus  that  we  can  obtain  a complete  extension  of  the  theory 

of  Section  9.  We  note  in  passing,  that  the  same  proof 

which  led  to  (11.13)  will  show  that  y^. 

h 1 

The  extension  of  the  Lax-Wendroff  scheme  will  be  just 

as  simple.  In  fact  (under  the  assumption  that  D,  6 and  f 
3. 

are  C ) the  difference  equation  for  (11.1)  will  have  the 
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same  linear  term  as  the  equation  for  (11.25),  (10.6) 
except  for  a perturbation 

(11.29)  k [D  + | ((B  - 6 )Do  + D (B  . - 5 ) + D2  + Dn) ] 

n 2 ch  n n n ch  n n 11 

as  one  can  easily  see  by  carrying  out  the  expansion  of 
Section  10  for  the  equation  (11.1).  It  follows  that  if 
k/hm  = 0(1)  (so  that  kllB^D  is  bounded),  Theorem  2 is 
applicable,  provided  k and  Cq  are  sufficiently  small.  One 
can  then  obtain  AP  sequence  solutions  and  carry  out  the 
convergence  argument  of  Section  10  (see  (10.17)  ff) , 
with  the  only  difficulty  being  the  verification  that 

(11.30)  II  (B2-B^)y#oo  = 0(hj) 

and  this  will  follow  immediately  if  one  assumes  a chain  as 
in  Assumption  I for  the  operator  L introduced  in  (10.16). 
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Appendix.  Proof  that  (L^,^)  Admissibility  Implies 
an  Exponential  Dichotomy. 

Here  we  would  like  to  give  a simplified  proof  of  the 
proposition  that  (A, A)  admissibility  for  the  inhomogeneous 
equation 

( A . 1 ) y = Ay  + f 

implies  an  exponential  dichotomy  for  the  homogeneous  equation 
(A.  2)  y = Ay 

It  is  assumed  of  course  that  A(t)  is  AP  although  we  will 
actually  do  the  case.  The  proof  is  valid  only  if  the 
underlying  space  E is  finite  dimensional  and  we  can  then 
regard  A as  a matrix  and  y and  f as  m vectors,  where  m is 
the  dimension  of  E.  The  reader  is  referred  to  Section  5 
where  an  entirely  similar  proof  for  the  difference  equation 
case  is  given. 

The  general  proof  is  given  in  M & S (p.  344,  Theorem 
103. A) , but  this  proof  relies  strongly  on  preceding  material. 
A finite  dimensional  proof  is  given  in  Coppel  (p.  134  ff) 
for  an  equation  on  [0,°°).  His  proof  can  be  extended  to  an 
equation  on  the  whole  real  axis,  but  the  proof  given  here 
is  simpler  and  more  in  keeping  with  the  ideas  of  Massera 
and  Schaffer. 

It  is  shown  in  M&S  (Theorem  103. A)  that  (A, A)  admissi- 
bility implies  ( ) admissibility  and  this  will  be  our 
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starting  point.  We  thus  assume  that  for  every  f in  Lm  , (A.l) 

has  a unique  solution  y in  Lm  and  the  estimate 

(A.  3)  lyl  < K«  ft 

will  hold  for  a certain  constant  K. 

Now  let  be  the  subspace  of  initial  data  which  gives 
rise  to  solutions  to  (A. 2)  that  are  bounded  for  i_  £ [0,°°). 

Let  S2  be  the  analogous  subspace  giving  rise  to  solutions 
that  are  bounded  for  t £ (-°°,0].  We  must  have  S ^ n S2=  {0} 
because  any  nontrivial,  bounded  solution  to  (A. 2)  would 
violate  the  uniqueness  requirement  of  our  definition  of 
admissibility.  Let  S 3 be  any  subspace  complementary  to 
© S2  so  that 

(A. 4)  E = Sx  © S2  © S3 

Let  P^,P2»P3  be  the  associated  projections,  i.e.  P^  is  the 
projection  onto  along  S2  © , etc.  We  have 

(A. 5)  I = Pl  + P2  + P3 

pipj  = 0 , i = 1,2,3;  j = 1,2,3;  i^j. 

We  will  first  show  S3  = { 0 } , i.e.  p3  = 0.  This  will 
follow  from  the  variation  of  constants  formula,  just  as  in 
the  difference  equation  case.  If  f has  compact  support  and 
y(t)  is  the  unique  L_r  solution  to  (A.l),  then  since  y is  a 
bounded  solution  to  (A. 2)  for  large  |t|,  we  must  have 
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(A. 6) 


Vo 


P3Y  (s)  f ( s ) ds 


Vo 


PnY  (s)  f(s)  ds 


Vo 


P-,Y  (s)  f ( s ) ds 


Here  Y (t)  is  the  fundamental  solution  to  (A. 2). 

Now  for  any  vector  z e e,  if  we  set  f = y Y ( t ) z 

L 0 f 1 ] 

then  ( A . 6b , d)  yield  immediately  P.z  = 0 and  we  can  conclude 


P^,  = 0 , i .e  . 


(A. 7) 


E = S1  © s2 


1 ~ P1  + P2 


Using  (A.6a,c)  we  also  see  that  for  any  f with  compact 


support  we  have 


(A.  8) 


y(t)  = G ( t , s ) f(s)  ds 


where 


(A. 9) 


G(t,s)  = 


Y(t)P1Y"1(s)  , 
-Y(t)P2Y"1(s)  , 


t > s 


s > t 


We  will  next  use  (A. 3)  to  show  that  for  every  fixed  t, 
G(t,s)  is  in  as  a function  of  s,  and  in  fact 
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(A. 10)  j ds  I G { t , s ) I £ K 

— 00 

where  K is  given  in  (A. 3) . This  is  a simple  consequence  of 
using  {A. 8)  as  an  operator  on  L^d),  where  I is  any  compact 
interval  to  conclude  that 


(A.  11) 


f ds  I G (t , s ) I < K 


and  then  letting  I ■*  (-00  , 00)  . 

Wa  now  will  show  that  (A. 10)  implies 


(A. 12) 


lG(t,s)  I =0(1)  . 


To  see  this  we  simply  observe  that  Y 1(s)  satisfies  the 
adjoint  equation 


(A. 13)  z'  * - Z A (s ) 

where  we  have  used  to  stand  for  d/ds.  Hence  fixing  t 

and  considering  the  region  s < t,  we  have 


(A. 14)  G(t,s)'  = - G(t,s)  A (s ) 


Now  since  A is  bounded  it  follows  that  G'  (t,s)  is  in  , as 
a function  of  s,  on  the  interval  (-°°,t]  with  an  norm 
which  can  be  bounded  uniformly  in  t.  This  of  course  means 
that  G(t,s)  must  approach  a limit  as  s -*■  and  since  G(t,s) 
is  in  (as  a function  of  s),  this  limit  can  only  be  zero. 
We  thus  have 


(A. 15) 


G(t,s) 


dr  G'  (t , r)  = 0 (1)  , s < t, 


— 00 
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and  a similar  argument  establishes  (A. 12)  for  s > t. 

We  must  now  show  that  (A. 12)  can  be  replaced  by  an 
exponentially  decaying  factor,  and  to  do  this  we  consider 
perturbations  of  (A. 1)  and  (A. 2);  namely  for  small  c > 0 
consider  the  equations 

(A. 16)  y = [A-eI]y  + f 

(A. 17)  y = [A-cI  ]y 

If  e is  sufficiently  small  we  will  have  (L  ,L  ) admissible 

00  00 

for  (A. 16) . This  is  a simple  consequence  of  the  contracting 
mapping  principle  analogous  to  the  argument  of  Section  4C  (see 
(4.*)  ff)  . The  contracting  mapping  principle  will  also  show  that 
the  admissibility  bound  for  (A. 16)  will  be  given  by  (see  A. 3) 

(A. 18)  K = — $ < K , |e|  < eft  , 

C (1-eK)  “ 0 

if  e0  '*'s  efficiently  small.  Thus  the  admissibility  bound 
for  (A. 16)  can  be  taken  independent  of  c,  if  e is  small 
enough . 

Now  the  fundamental  solution  to  (A. 18)  is 

(A. 19)  YC  (t)  = e“Ct  Y (t) 

. £ £ 

and  if  P^,  P2  denote  the  projections  replacing  P^  and  P2  , 

£ 

while  G ( t , s ) is  the  new  Green's  function,  then  the  analysis 
which  led  to  (A. 12),  and  which  depended  only  on  admissibility, 
is  equally  valid  here  a*  d we  can  conclude 


-108- 


(A. 20) 


lGe(t,s)l  = 0(1)  ; lGe(t,*)l1  = 0(1)  , 

where  the  bound  in  (A. 20)  can  be  taken  independent  of  e 
(because  of  (A . 18) ) . 

The  estimate  (A. 20)  means  in  particular  that  for  small 
e there  is  a constant  K 2 independent  of  e such  that  if  s>  t, 

(A. 21)  I Y (t)P2Y_1(s)  I < K2  e“e(s“t) 


and  it  is  only  necessary  to  show  that 

pi  - pi 

(A. 2 2)  1 1 


To  prove  (A.  22)  we  let  and  S2  replace  S.^  and  S2 
It  then  follows  from  (A. 19)  and  e > 0 


(A. 23) 


si£si 
S2  £ S2 


and  (A. 22)  will  be  established  if  we  can  show  that  the  ranks 
of  the  stable  and  unstable  projections  are  unchanged  for 
small  e,  and  this,  in  turn,  will  follow  from 


(A. 24)  «P^  - PLI  = 0(e)  . 

To  obtain  (A. 22)  we  set  He(t,s)  = Ge(t,s)  - G(t,s). 
Now  for  any  f with  compact  support  we  let  y and  y£  be  the 
unique  bounded  solutions  to  (A.l)  and  (A. 16)  respectively. 
If  w£  = yf - y we  see  that  w£  is  the  unique  bounded 
solution  to 

(A. 25)  w£  = Aw£  -e  y£  , 
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and  we  can  conclude  from  (A. 3)  together  with  its  analogue 
for  (A. 16)  that 


(A.  26) 


I w I = 0(e  ) II  f II 

E OD  ' 00 


Wow  w.  can  also  be  written  as 

00 


(A. 27) 


w£(t)  = 


r 


j 

“00 


ds  (t,s)  d (s ) 


and  using  the  same  argument  that  led  to  (A. 10)  we  see  that 

00 

t 


(A. 28) 


J 

— 00 


ds  I H (t,s) • = 0(e)  . 

(Note  that  the  bound  in  (A. 28)  will  be  independent  of  t, 
although  we  will  only  use  it  for  t = 0.) 

Now  if  s < t we  see  that  as  a function  of  s,  He(t,s) 
satisfies  the  equation 

(A. 29) 


He  = - HeA  + e Ge (t,s)  . 


This  equation,  together  with  (A. 20)  and  (A. 28),  shows  that 


(A.  30) 


ds  II  He  ( t , s ) ! = 0(e) 


and  since  H (t,s)  is  in  as  a function  of  s we  can  write  for 
s _<  t 
(A. 31) 


HL(t,s)  = I ds  H£  ( t , s ) = 0(e) 


and  if  we  set  t - 0 , s - 0_  , (A. 31)  yields  (A. 24),  which 
proves  one  part  of  the  requirement  of  an  exponential  dichotomy 
and  the  other  part  follows  on  replacing  -e  by  +e. 
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